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Linearized Treatment of Supersonic Flow 


Through Ducts" 


H. F. LUDLOFF# ano F. REICHE** 


New York University 


ABSTRACT 


The object of the investigation is to determine the velocity 
and pressure field occurring in supersonic flow through ducts of 
various shapes. In order to obtain results in analytical form, 
giving approximate insight into essential flow features, linearized 
treatment has been applied. 

The solution has been composed by assuming that 
sonic sources’’ of various intensities are distributed over the sur- 
face of a skeleton cylinder. In Part (II) the expressions for ve- 
locity and stream function have been derived for three cases: 
(A) for constant source intensity; (B) for source intensity that 
increases linearly along the generatrices of the cylinder; (C) for 
a certain combination of the two previous cases, which turns out 
to be the most satisfactory approach to the problem insofar as 
generality and validity of the procedure is concerned. This last 
approach yields nozzle shapes for which location and size of the 
throat can be adjusted at will. Moreover, a general theorem 
has been derived enabling one to give expressions for the velocity 
field of any arbitrary source distribution (Appendix IT). 

Particular attention has been paid to the physical discussion 
of the results (Part III). The resulting pressure distribution 
inside the duct points to the occurrence of compression and ex- 
pansion waves that appear at the surface of the Entrance Mach 
Cone and at the extension of this surface, called the ‘‘second 
Mach cone.”’ In general, the intensity of these waves is small 
and the distance between them appreciable. The computations 
indicate, however, that, with certain properties of the duct rim, 
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a compression wave would be immediately followed by an ex- 
pansion wave, while the wave intensities increase. 

The form of the wave front has been computed in a second 
approximation. For small lip angles, the wave shape is only 
slightly different from that of a double cone, in qualitative agree- 
ment with Ferri’s analysis. 


NOTATION 


= axial and radial coordinate of test point P 
axial coordinate ‘of source element 
radial distance from test point P to source 
element 
radius of skeleton cylinder 
azimuth position about duct axis 
Mach angle; M = free-stream Mach Num- 
ber; 8 = tana 
= lip angle of duct 
inclination of shock front 
distance between sources and sinks in case](C) : 
see Fig. 3 
x/Bd; Y =r/d; ¢ = §/Bd; 
Z-—(i-— Y) 
sd 2Y 
= (1+ Y)?; b = (1 — Y)*;-¢ = Z? 
general source intensity per unit area [L/T] 
intensity of constant sources [L/T] 
intensity of linear sources [L/T] 
disturbance velocity potential [L?/T] 
disturbance stream function [M/T] 
disturbance density [/L*] 
free-stream potential, stream function, and 
density 
U = free-stream velocity 
vz and v; axial and radial disturbance velocities 
o ABr/U ~ 6/a; o' = A’Br/U 
p absolute pressure 
Cp = pressure coefficient 


; Z' = (x — B)/pd 


®o, Vo, and po 


(I) INTRODUCTION: OUTLINE OF THE INVESTIGATION 


I VARIOUS AERODYNAMICAL PROBLEMS, it is desirable 
to have an approximate, analytic method for com- 
puting the supersonic field of flow through ducted 
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bodies in order to get quick insight into essential flow 
features as related to the shape of the nozzle wall. 
While the two-dimensional flow through nozzles has 
been commonly treated by the “‘method of character- 
istics,’ it has appeared that, for the three-dimensional 
case, the cumulative error (which must necessarily arise 
in this step-by-step method) makes it worth while to 
try a different approach. The advantage of the present 
method is that it would yield the streamlines all at once 
and at the same time permit some insight to be gained 
of the analytical features of the solution. 

We are considering a supersonic air stream entering 
a nozzle whose constriction occurs slowly enough so 
that the prerequisites of the linearized theory can be 
considered as fulfilled. The flow may be anticipated 
to remain supersonic throughout, in spite of the occur- 
rence of one or two oblique shocks. The approach to 
the problem is analogous to that in the problem of 
supersonic flow about slender bodies of revolution, 
treated by von Karman and Moore.! In the case 
at hand “supersonic sources” of various intensities are 
placed on the surface of a skeleton cylinder, and the 
velocities, pressures, and streamlines resulting from 
such a ‘‘source’’ system are determined. According to 
the usual procedure, an appropriate system of stream- 
lines may be considered to form the wall contour of the 
nozzle. Therefore, the present objective is to determine 
the wall contour and the velocity and pressure field cor- 
responding to a given “‘source’”’ distribution. The in- 
verse problem will be treated in a second paper. 

The influence of the nozzle wall on the velocity field 
can be best understood by dividing the whole flow field 
into four districts (see Fig. 6). Starting at the duct 
entrance, there is first a conical region (1), in which the 
air stream continues to flow undisturbed. Behind the 
Mach cone, originating at the entrance rim, the disturb- 
ances arising from the duct wall are felt. But they are 
not felt all at once. In the two intermediate districts, 
(2) and (3), the influence of the upper or the lower 
part of the wall section will predominate. From Fig. 6, 
one can see that.a point in districts (2) or (3) is not in- 
fluenced by source elements of the directly opposite 
part of the nozzle wall. Disturbances, produced by the 
entrance rim, are propagated along Mach lines whose 
extensions form a second Mach cone. Behind this 
second Mack cone begins the rear district (4), in which 
any point of the flow field is under the combined influ- 
ence of source elements located at any value of the 
azimuth w, on the surface of the cylinder. 

The analytical investigation is divided into three 
parts. In Part (A) the source intensity is assumed to 
be constant all over the surface of the skeleton cylin- 
der. Asa second case, in Part (B) the source intensity 
is assumed to increase linearly along the generatrices 
of the cylinder. Finally, Part (C) contains the analysis 
of a source distribution, which is a certain combination 
of the two previous cases and, apparently, the most 
satisfactory approach to the problem, insofar as flex- 
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Fic. 1. Wall contours and streamlines of nozzle produced by 
constant source distribution. 


ibility, generality, and validity of the procedure is con- 
cerned, 

In the case of constant source intensity, it turns out 
that the streamlines and wall contour of the nozzle form 
a throat at some distance behind the entrance rim, and 
the inside diameter increases at the exit section to the 
same magnitude that it had at the entrance section. 
Since velocity, density, and pressure are also the same 
at the rear end as at the entrance, the wave drag origi- 
nating from the inside of the duct will be zero. 

Some quantitative results of the analysis are shown in 
Fig. 1. The velocities and streamlines of the flow pro- 
duced by sources of constant strength have been com- 
puted and three particular streamlines that start at the 
rim and represent possible wall contours have been 
plotted. The wall contours correspond to three differ- 
ent source strengths as indicated in the figure, and, of 
course, the stronger the source intensity, the larger the 
deflection of the wall contour from the skeleton cylinder 
and the more pronounced the throat. The stream- 
lines entering halfway between axis and duct wall are 
also shown. 

The analysis permits one to draw certain conclusions 
regarding the occurrence of shocks. At the entrance 
cone a discontinuity of velocity appears which is small 
except at the vertex of the cone. [See Eqs. (11) and 
(39).] In other words, the Mach cone at the entrance 
of the duct represents a conical compression shock that is 
weak except in the intermediate vicinity of the vertex. 
It is possible to compute the form of the shock front 
in second approximation by use of the shock polar 
equations. The result obtained in this way can be 
compared with Ferri’s computation of the entrance 
shock carried out by means of the method of character- 
istics and is shown to be in satisfactory agreement. 
[See Section IIT, (3).] 

Furthermore, it is found that, on the surface of the 
second Mach cone, a compression and subsequent ex- 
pansion wave originate, one immediately behind the 
other. This result is discussed in Section III, 2, and 
it is pointed out that a similar phenomenon may occur 
in the problem of von Karman and Moore. Some 
features of this result are outside the range of validity 
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of the linearized theory—-which is to be expected, since 
at the entrance rim of the cylinder the source intensity 
increases discontinuously from zero to a finite value. 
But it can be shown that this pattern constitutes the 
limiting case of a more general flow pattern which is 
perfectly within the range of this theory. Therefore 
our results in the vicinity of the second Mach cone are 
qualitatively correct, although they may not be inter- 
preted quantitatively. 

In accordance with the theory of characteristics, the 
singularity occurring on the second Mach cone is due 
to properties of the duct rim. Therefore, the features 
of an air stream entering a duct with a knife-sharp rim 
have been analyzed. (See Section II, B.) Flow con- 
ditions of this character can be materialized if the sur- 
face of the skeleton cylinder is covered by a source dis- 
tribution whose intensity is linearly increasing, starting 
at the rim. And indeed, the analytical treatment shows 
that this velocity pattern does not produce any pecu- 
liarity at the second Mach cone. There are also no 
discontinuities of velocity at the entrance cone in this 
case, so that the whole flow field can be considered to 
be within the range of validity of the linearized theory. 
Disregarding boundary-layer effects, the sharp entrance 
rim of a duct would seem to be a means of avoiding 
shocks, while a blunt rim will produce them. 

Computation of the streamlines that represent the 
contour of the nozzle wall (see Fig. 2) shows that in 
this case no throat is produced, but the cross section 
of the duct narrows down the more one recedes toward 
the rear. Since the cross section decreases in a loga- 
rithmic way, the flow becomes, for practical purposes, 
parallel after a certain distance from the entrance is 
reached. However, the ratio of exit to entrance section 
will be definitely smaller than one, so that a consider- 
able wave drag would result. 

It is possible to combine the merits of the last two 
nozzle shapes without encountering their respective 
setbacks. Fig. 3 shows the resuit of this combination 
which consists of a source intensity that starts to in- 
crease linearly at the rim and remains constant after 
a certain distance B from the entrance is reached. 
Careful analytical investigation (see Section II, C) 
shows that the knife-edged rim hinders the formation of 
a discontinuity of velocity at the entrance or a singu- 
larity at the second Mach cone. On the other hand, 
the constant part of the source layer takes care that the 
inside diameter, after narrowing down to a throat, 
increases at the exit section to its original magnitude 
so that the wave drag is zero. The source distribution 
that produces this nozzle shape can be resolved into 
two linearly increasing layers, one consisting of sources 
and starting at the entrance rim of the duct and the 
other consisting of sinks and starting a distance B be- 
hind—i.e., where, in the previous representation, the 
constant layer begins (see Fig. 3). As linearly increas- 
ing source intensities have been shown to produce flow 
fields that are everywhere regular, the super position 
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Fic. 2. Wall contours of ducts produced by linearly increasing 
source distribution. 
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° 5 ' Ss 2 25 3 35 4 


Z= “/ap 


Fic. 4. Wall contours of nozzles produced by two different 
source sink combinations. 


of two such distributions will certainly not produce any 
singularity in the flow field. By varying the source 
and sink intensities and by changing the distance B be- 
tween sources and sinks, the location and size of the 
throat can be adjusted at will. (See, e.g., the two wall 
contours of Fig. 4.) Thus, the use of this source dis- 
tribution (C) represents a flexible procedure that pro- 
vides results that are valid at all points inside the duct. 
The first two cases (A) and (B), may be considered as 
instructive auxiliary considerations. 
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(II) ANALYTICAL INVESTIGATION OF THE FLOW FIELD 


(1) General Form of the Velocity Potential 


The analytical investigation is based on the compu- 
tation of the velocity potential, from which the velocity 
components, the stream function, the pressure, and 
the wave drag can be derived. Consider ‘‘sources’’ 
of given intensity per unit area distributed over a 
circular cylinder upon which a supersonic air stream 
is impinging. The resulting flow pattern will consist 
of the undisturbed air stream plus the perturbation 
flow produced by these supersonic ‘‘sources.’’ Of 
course, since one is dealing with a linearized theory, the 
perturbation velocities are assumed to be small in com- 
parison with the main flow velocity; this is the case if 
the cross section of the duct changes slowly. The per- 
turbation flow can be completely described when the 
velocity potential due to this source distribution is 
known. 

As can be seen from Fig. 5, the potential, d@ at point 


P, contributed by the elemental source area at x = &, 
is given by: 
£) dg dw 
db = fC (1) 


“i: — £)? — BR? 


here f() is the source strength per unit area at x = &; 


the source distribution is axial symmetrical, so that 
f does not depend on w. 8? = M? — 1, M being the 
main stream Mach Number. R is the distance from 
the point P to a surface point in the same plane, as ex- 
plained in Figs. 5 and 7. The denominator.in Eq. (1) 
is a “pseudo distance,’’ which reduces to the real dis- 
tance between P and the surface element when ./ = 0. 

To obtain the total potential at P, an appropriate 
integration over the surface of the cylinder must be 
made. It is apparent that only sources contained be- 
tween the entrance rim and the fore cone of P will af- 
fect the potential at P (see Fig. 5). Therefore, the in- 
tegration over £ extends from 0 to x — BR, where x is 
the abscissa of P. 




















Fic. 5. Geometry of skeleton cylinder and fore cones. 

















Fic. 6. The four districts of the flow field: (1) zone of silence; 
(2) and (8) front district; (4) rear district. 








Fic. 7. View of entrance plane of duct, showing intersection of 
skeleton cylinder and Mach cone of P. 


To determine the limits of w, two cases may arise 
(see Fig. 6). 

When point P is in region (4), called the 
trict,” it is obvious that w runs from —7rto +7. How- 
ever, when P is in regions (2) or (3), called ‘‘front dis- 
trict,’’ the fore cone of P intersects only a portion of 
the cylindrical surface so that w has the limits —ws and 
+wx. (See also Fig. 7.) 

With the understanding that wx = m when P is in the 


“rear dis- 


rear district, the expression for the total potential at P 


becomes: 
-a f” x — BR f(g) 
®@ = —2d ‘iw ff dé : 
"Wis = _ ¢)? — p? R? 


From Fig. 7 it can be seen that 


R? = d? + r? — 2dr cosw (3) 





D: 
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and from Fig. 6 it follows that for w = ws the é-integra- 
tion interval reduces to the point = 0. Hence, 


0= &=x — BR(ws), R(wx) = x/B (8a) 


Thus, using Eq. (3), 
wx = cos! [d? + r? — (x?/B*)]/2dr (4) 

By differentiating Eq. (2) with regard to x and 7, one 
can show that the linearized equation determining an 
axial symmetrical potential ® is fulfilled. While the 
direct verification is lengthy and cumbersome, the proof 
given in Appendix I is more elucidating. 

We have now to carry out the integration in Eq. (2); 
therefore the treatment will be divided, according to the 
type of source distribution considered. 


(2) Case A. Constant Source Distribution on Circular 
Cylinder 
Substituting in Eq. (2) f(£) = A and carrying out the 
t-integration, one gets: 
a [a ! nn (5) 
@ = —ZAC w 1lo ee > 
0 81 — V1 — (6R/x)? 


Differentiating with respect to 7, one gets the radial 


velocity component: 


: fi r — dcosw 
v, = 2Ad dw —_—— 1 Te 
J 0 R? V1l- (BR/x 


It is appropriate now to introduce the dimensionless 
variables: Z = x/Bdand Y = r/d. 
Upon introduction of a new integration variable, 


p = R?/d? = 1+ FY? 


— 2V cosw 
Eq. (6) becomes: 


v9 = AZY-! [ dp t bss Vab 
a6 sbi ‘op Ve — p)(a — p)(p — B) 


« 


(7a) 


forb <c <a: front district; 


/ 
"@ di — Vab 
v, = AZY- y| . _ = === (7b) 
b p Vc — p)(a — p)(p — bd) 


for) < a < ec: rear district. 
Here, a = (1 + Y)?,6= (1 — Y)*,andec = 2?, 
If point P moves on the second Mach cone, ¢ be- 
At the integration limits, the inte- 





comes equal to a. 
grand in Eq. (7a) or (7b) becomes infinite as 1/V p —b 
(or 1 IV ¢ — p and 1/Va — p, respectively) ; this square 
root singularity of the integrand disappears in the 
integrated expression for v,, However, if c — a, the 
two square root singularities of the integrand coalesce 
to form an 1/(a — p) singularity. Then, the integrated 
expression for v, becomes logarithmically infinite. 

(a) Front District—In order to obtain a constant 
integration interval that is independent of Z and Y, one 


may introduce: 


p=b+(c—5)n 


sothatO < y»<1,whendB<pc<ec. Putting 


2 2 =~ ar 
a-b 4Y 








2 = 


which is smaller than 1 in the front district, Eq. (7a) 
becomes: 


1 
uv. = -azy-* [ar x 
0 


_ [((1 — Y) — 2k*y] 

[(1 — ¥)? + 4¥&%]Vn(1 — n)(1 — Bn) 
It can be shown that Eq. (8) can be expressed as a 
combination of elliptic integrals. The transformations 
required for that purpose will be explained in a paper by 


Mack and the authors.? Here, however, the solution 
will be carried out on the basis of an expansion of 1 + 


Wi ~ k*n, giving: 


r—1 1 (1 — Y) — 2k’ \ 
v, = —AZY~"* d { 2 ei 
; J, "la — yy? + 4VEnf * 
1 + (1/2)nk? + (3/8)n°k* +... “a 
Vn(1 — 7) 


Combining equal powers of » under the same in- 
tegrals, one gets: 


v, = -AZY~“* [11 — Yh — 
(1/2)(38 + Y)R72, — (1/8)(5 + 3VY) Rh — 
(1/16)(7 + 5Y)k*2, — (5/128) X 
(9+ 7Y)RI, —...] (10) 





(8) 











where 
nm — (1/3) 
n 


1 
ey Xe hi 
Jo (1—n)“[(1 — VY)? + 4k Yq] 





(10a) 


which can be evaluated by means of a recursive formula 
and yield simple rational functions of Y, Z, and k?. 
Upon substituting J, and k®", Eq. (10) becomes: 


v, = —voa-ty-*f1 — o (SEY) 
4 
uh + SY 2 + 2(1 — 2] i 
a/2) ¢ (2+ ; 


(9/32)s8 (ee) x 


E + 3Z(1 — +t (8/3)(1 — em) - \ (11) 








Here two more characteristic dimensionless parameters 
have been introduced—namely, o = AfSmr/U ands = 
[Z — (1 — Y)]/2Y; the lines, s = constant, represent 
the set of rays going through the tip of the Mach cones. 
The meaning of ¢ is explained later. 

Stokes’ three-dimensional representation of the 
stream function y, which can be derived by considering 
the flux through a cylindrical surface,* may be written 
as either 
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° * Zz 
vi=i-m f rv, dx’ = —mver f dZ’' X 
B(d — r) 1-¥Y 


(8U-) Yv,(Z’, Y) (12a) 


or 
v= —po6? fA r'v, dr’ = —poUd? X 
0 
Y 
f dY'(8°U>) Y',(Y,’Z) (12b) 
0 


By substituting Eq. (11) into Eq. (12a) and carrying 
out the integration, one gets: 


1\ 3 Y 
yy = WioY*/.s {1 — (5)s2+ $- ee x 








5+ 3Y 
 eeamidiae ) sae 
Ss 
15 7+ 5 1 
—(1— Y — —... : 
16 ( r 12 r. ; ‘ (13) 
here, Wo = (1/2)poUd? is the flux crossing the entrance 


section of the nozzle. 

The radius of convergence of the series (13) is s = 1, 
since v,, from which Eq. (13) has been obtained by in- 
tegration, becomes logarithmically infinite for this 
value. It can be shown that up to s = 0.8 four terms 
of Eq. (13) provide, for practical purposes, sufficiently 
accurate results, as may be seen from comparison of 
Eq. (13) with the exact solution in terms of elliptic in- 
tegrals (see references 2 and 4), which yields: 


Y = W210 {44/Y [k’ 2K(k) — E(k)] + 
4+Z[e + 2(K(k) — E(k)) F(R’, 0) — 
2K (k)E(k’, 8)]} (14) 
here: k’? = 1 — k?; 6 = sin=' [(1 — Y)/Z] (0 < 
6< 2/2). 

Before passing over to the determination of the 
stream function in the ‘‘rear district,’’ it is of interest 
to examine the stream function on the dividing line s = 
1—i.e., on the second Mach cone. This enables one to 
obtain a value for y at a point where both the front and 
rear representations fail. As explained later on, the 
exact value of the stream function on the second Mach 


cone assumes the form: 


(1 + Y) are tan (~/Y) — VY] 
(15) 


The closed form expression (15) for Y becomes possible, 
because, for s = 1, the previously mentioned elliptic 
integrals degenerate. Numerical evaluation of Eq. (13) 
for s = 3/4 and Eq. (15) for s = 1 has shown that the 
respective values are compatible. 

(b) Rear District—Analogously, as in the front dis- 
trict, one may now introduce 


= b+ (a—bd)n 


Vs oi = SWoonr- 


and 
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i. «$s 4Y 
~ Z24 (1 — V)? 


which is smaller than 1 in the rear district. Then Eq. 


(7b) yields: 
1 
0 


[1 — Y) — 2n][1 + (1/2)nk,? + (8/8) n?kit + . 


Vx(1 — — m[(l — Y)? + 4¥n] 


v, = —A2Z[Z? — (1 — Y)?] 








(16) 
which can be evaluated similarly as Eq. (9) and is shown 
to vanish asymptotically for large Z as Z~?. 

To obtain the stream function, one may use Eq. 
(12a) in a form more suitable for evaluation in the rear 
district. Since 


Vin0 = — mf rv, dx’ — » fi 1U; ax’ = 0 
B(d —r) tv 


one has 
v(x, 7) = -m fom, dx’ (17) 


Substituting Eq. (16) into Eq. (17), one obtains: 
- 12 + OY — 3¥° 
= Wolk — — = 
v (1/3)hs 16(1 + Y) 
” 75 + 45Y — 25Y? + 5Y3 t 
(1/5)&,*| —— - - 
/syks| ie |+... (18) 
In terms of elliptic integrals (see references 2 and 4), one 


obtains: 
4/fVii-—-Y|, ,-< . 
y = W2r 7) v k + yt *K (ki) — E(k | -++- 
— E(k;)) F(ki', 0) — 2K (ki) E(hi’, ait 








ky 
Z[m + 2(K(hi) 
(19) 


here: 6 = sin! ((1 — Y)/(1+ Y)]; &’? =1 — R,?. 

(c) Rear District: Expanston After p/Z*.—A more 
convenient representation of the stream function in the 
rear district is obtained by starting from expression (6) 


for v;: 


v, = aya fr dw p-(U- —( - ===) 
0 p vi - p/Z? 


expanding 1/V1 -~ p/Z?, in powers of p/Z?, one obtains: 


Ls = 73 
v7, =A yf dw E 2 a x 
0 p 
1 ae 3 OR 
[1 + (5)eo/z + (5)ie/z2) +.. | (20) 


Evaluation of the integrals f dw p“(w, Y) yields simple 
0 


polynomials of 2mth degree in Y. 
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To obtain the stream function, one may use again 
Eq. (12a), which yields: 


y = WH2V*Z[1 + (1/4)Z-7°(2 + Y?) + 
(1/8)Z-4(3 + 6Y? + Y*) + (5/64)Z-* X 
(4 + 18Y? + 12Y*+ Y*%+...] (21) 


this becomes identical with Eq. (18) for large Z. 

The stream function of our flow pattern is obtained 
by adding the undisturbed and perturbation stream 
functions at any ~~ (Z, Y): 


Vv = 24 u(Z, Y) = W(d,/d)? (22) 


The left-hand side of this equation represents the total 
stream function; the right-hand side represents the 
constant value that the stream function assumes on the 
particular streamline which, in front of the duct, is at 
a distance d, from the axis. To plot the streamlines, 
Eq. (22) has to be solved for Y in terms of Z. The 
streamline representing the wall contour is obtained by 


putting in Eq. (22) d, = d. This yields: 
Y? = 1 — (1/W)¥(Z, Y) (22a) 


It may be noted that the lip angle 6, which this wall 
contour forms with the skeleton cylinder, is given 
by: 


O)s ao = 1A /U (23) 


y=1 


6 ~ tan 6 = (v, 


The dimensionless quantity = 7A6/U, introduced 
earlier turns out to be 


rA 1 _ tan a 6 


U tana 





(24) 


~ tan a a 


and measures the lip angle in terms of the Mach 
angle. 


3) Case B: Source Distribution of Linearly Increasing 
Strength 


In the following: 
f(§) = A’t/Bd = 


here: ¢ = £/@d. 

A’ isa constant which characterizes the rate of source 
intensity increase. Substituting in Eq. (2), one gets, 
for the velocity potential: 


—BR yr i 
> = -24 f~ dw i: dé Vx = ae BR (25) 
3 _— - 


Introducing the dimensionless variables, as before, and 
in addition, 





= A’Br/u; Ud = 


one can write: 


= —yr—!20'B- fm ao fa e ~pnilad = 
; 'Viz—p 


which, on carrying out the ¢ = integration, yields: 
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© = —Gr- a-20' f dw X 
0 


| 2 os Z+VZ—>p 


p 


- VE - q (26) 


Differentiating with regard to 7, one now obtains: 


we . Jo 
= Uris-20'f f dw | vz" = | + 
( 


dw 
>| aes (27) 


As shown previously, p(wx) = Z’, so that the last of 
the terms in Eq. (27) vanishes. 

Again using p instead of w as integration variable, 
Eq. (27) transforms into: 


v. = Uripw'ys f dp Ve = a = oo = was _Vab) (28a) 
> P Va (a — p)(p — b) 


front dintrlet: 


z 
| -2 cosh! my + Vz? = 
p 


forb<c<a: 


Ux-—B ays ["? Ve - ele — Vab) (28b) 
h 


p V(a — p)(p — b) 





v, = 


forb <a <c: rear district. 
In marked difference to the case of constant source 
distribution, Eqs. (28a) and (28b) contain the factor 


Ve = p in the numerator instead of the denominator, 
so that in the limit c > a@ the factors Vc — p and 








Va — p cancel rather than form a cumulative singu- 
larity. 

(a) Front District—Proceeding similarly as in the 
case of constant sources, one obtains without difficulty 
series expansions for all the quantities under study. 
However, by means of the Recurrence Theorem, formu- 
lated in Appendix II, the derivations may be consider- 
ably shortened. It is proved in Appendix IT that 


o 
9!" = OFPiin = (dg)— (=) ren (A20) 
Ox o 
Viin ia “f., dL" conse, (A21) 


Hence, using Eqs. (A19) and (13), one obtains im- 
mediately: 

; l 3+ Y 

ie.  690— ry'/e, _ - a= 
v7, = —2B-'U0'!) {1 (3)s 4 


5+ 3Y 1 . 
(1) BA8F Gy (1) 
[v2 43" (2) a= vy EEF], . (29) 


Since Eq. (29) is convergent for s < 1, one may integrate 
termwise, so that according to Eq. (12a): 
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(b) Solution on the Second Mach Cone.—From the 
general expressions [Eqs. (28a) and (28b)], one can see 
that on the second Mach cone (a = c) the integral 
simplifies to: 


(»,) Ur-1p-1 cad Mek: —— 
)y Sean = Tv Co - — ——————— 
tao b Pp Vp —b 


which can be evaluated in closed form—namely, 


1 Y l 
(v,)x2ay = —U4r-'!8—0' (= are tan 1/Y- <a) 
(31) 


It is not possible to derive a closed expression for yj, at 
the second Mach cone; in order to make use of Eq. 
(12a), one needs v, as function of x, which is not avail- 
able. However, a simple expression for y,,,,,. at the 
second Mach cone can be derived from the above by 
direct use of Eq. (A19). This expression has already 
been quoted in Eq. (15). 

(c) Rear District—To compute the stream function 
in this case one may apply the recurrence theorem 
(A21) and obtain: 


. i\f1. ; 
Viin = Wo20'} og 2Z — (5) E (2 + } 2) +4. 


[4 


WA wcguig wea (2\4 
(G)pG+o + ry + (2) x 


(4 + 18Y? + 12Y4 + Y*) +.. | (32) 


As in the case of constant sources, so here the stream- 
lines can also be evaluated by using the relation (22a): 


Y? = 1 — (1/Wo)Wjin(Z, VY) 


By substituting the stream functions (13), (21), (30), 
and (32) into (22a), it becomes possible to compute 
shape and location of streamlines and wall contours. 
The results of such computations are shown in Figs. 1 


and 2. 


(4) Case C: Sources Distributed with Combined Linearly 
Increasing and Constant Strength 


In order to adapt the method to the representation of 
flow through more general nozzle shapes, an instructive 
combination of the source distributions treated in Parts 
(A) and (B) shall be considered in this section. This 
type of source combination is capable of generalization; 
at the same time it permits one to retain those par- 
ticular features of constant or of linear source distribu- 


tions which are desirable while on the other hand 
eliminating the undesirable ones. 

In Fig. 3 the shaded heavy line represents the con- 
sidered distribution of source strength along a generatrix 
of the cylinder: the source strength is linearly increas- 
ing between £ = O and £ = B (B being some arbitrary, 
fixed distance behind the rim), and then it remains at 
constant strength from § = Bto —§~ ~. There is 
again the same source strength for all azimuth positions 
for any particular value of &. 

In order to form the velocity potential of this source 
combination, it is, of course, possible to compose it, 
as indicated—namely, by inserting into Eq. (2), first, 
f(é) = A’(&/Bd) and extending the £-integration from 
£ = Oto & = B; second, f(t) = A’(B/8d) and extend- 
ing the £-integration from § = B to & = x — BR; fin- 
ally, by adding the two expressions. It can be shown 
that this procedure of composing the potential yields 
in either component potential a singularity; however, 
in the combined potential the singularities cancel each 
other; this indicates that, for this source combination 
(also when the point under consideration moves on the 
second Mach cone) the velocity field will remain regular. 

These features become evident in a more convenient 
way if the prescribed source combination is formed in 
another manner, which is also indicated in Fig. 3 by 
(a) having sources of strength A’(é/8d) start at & = 0 
and extend indefinitely to the right and (b) having sinks 
of the same strength start at £ = B and extend indefi- 
nitely to the right. Then the net source strength for 
& > Bis given by: 


f(&) = A’(&/Bd) — A’(é — B)/Bd = A’'B/Bd 


as it should be. 

In this way one obtains the velocity potential by 
algebraic addition of two analogous expressions for 
linear sources and linear sinks, with the following pro- 
visos: replace in Eq. (25) A’ of the sources with (— A’) 
for the sinks; also replace x of the test point and é of the 
source element, corresponding to the sources, with x’ 
and ¢’ for the sinks, where 


x’ =x — B; ’=:-—B (33) 


This stipulation expresses the fact that we are, in ef- 
fect, dealing with two cylinders: one on which the 
linear sources are located, and another one that is 
covered with linear sinks, the two cylinders being co- 
axial and identical except that their entrance rims are 
shifted with regard to each other by a distance B. 
With this interpretation in mind (see Fig. 3), one can 
visualize that the source distribution begins at & = 0, 
w extending from —w»x to +ws, and that the sink dis- 
tribution starts at t’ = O—i.e., § = B, w extending 
from —ws’ to +wx’. Here wx, defined by Eq. (4), is 
the azimuth of the point at which the Mach circle and 
the rim of the source cylinder intersects; wx’, corre- 
sponding to the point of intersection between Mach 
circle and rim of sink cylinder, is obtained by replacing 
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in Eq. (4) x by x’. 
potential of this source sink combination results: 


@ = b+ + o- 


we o—be A’(§/Bd) 
» <n a 1§ SS 
af is f “OVeo bp —BR 


fa i, fit pee ee 
2d fo dw ag V (x! — £2 si B°R? (34) 


where 


ma = oot 


, Sf pe = (34a) 


wx’ = cos 
Z’ = (x — B)/Bd = Z— Z| 


As indicated by Eq. (34), ®, v,, and y for the source 
combination can be obtained by superposing two anal- 
ogous expressions for linear sources, as previously de- 
rived. 

Since either one of these expressions is regular in the 
whole flow field, the combination is also regular. 
Moreover, the combination can be shown to remain 
finite for Z ~ ~, as pointed out at the end of this sec- 
tion. 

If one wants to use explicit expressions for Y, one has 
to distinguish between the arrangement in which the 
sinks start in front of the second Mach cone of the 
sources and the arrangement in which the sinks start 
behind. (See Figs. 3 and 4.) 

As an example, only the first arrangement shall be 
discussed here. Depending on the various possible 
positions of the test point P, one may distinguish six 
cases (see Fig. 8): 














Fic. 8. Districts of flow field occurring in Case C. 


Case (1): 
x < B(d — 1) 


Then: ®= 0. 
Case (2): 


\8(d + 1r), if 28r < B 


B(d — r) <*<13 + pd -—v), 26> B 


Then ® is determined on the basis of source elements 
alone, and these corresponding to the ‘front district” 
only. 


Therefore, the following velocity 


Case (3): 
B+ Bid —1r) <x < B(id+r) 


where 287 > B. 

Then ® is determined on the basis of source and sink 
elements corresponding to the respective front dis- 
tricts. 

Case (4): 


Bid +r) <x< B+ Bd —1r) 


where 267 < B. 
Then ® is determined on the basis of source elements 
alone, but now these correspond to the rear district. 


Case (5): 


B+ B(d — 71), 
if 28r < B 
B(d + 1r), 

if 28r > B 


<x < B+ Ad+7) 


Then ® is determined on the basis of source elements 
corresponding to the rear district and also on the basis 
of sink elements corresponding to the front district. 

Case (6): 


B+ B(d+r)<x 


Then ® is determined on the basis of both source 
and sink elements corresponding to their respective 
rear districts. 

Furthermore, as previously explained, ws = 7, if 
x > B(d +r), and ws’ = 2,ifx > B+ Bd+7r). 

Thus, Case (1): y = 0; Case (2): yw given by Eq. 
(30); Case (3): yw given by a difference of two repre- 
sentations, analogous to Eq. (30), replacing in the 
second one s = [Z — (1 — Y)]/2Y by s’ = [2Z’ — 
(1 — Y)]/2Y; Case (4): W given by Eq. (32); Case 
(5): y obtained by combining Eq. (32) and a represen- 
tation derived from Eq. (30) by replacing Z by 2’; 
Case (6): ¥ given by a difference of two representations 
analogous to Eq. (32), replacing in the second one Z by 
Z’. It is important to note that by taking the differ- 
ence, the logarithmic infinity of y for Z — ©, as indi- 
cated in Eq. (32), disappears; for extremely large Z- 
values, one may write: 


We = Wo20’¥*[(Zo/Z) +... ] (35) 


here Z) = B/ Bd. 

Using the representations of y in districts 2, 4, 5, and 
6, two different wall contours have been computed and 
plotted in Fig. 4. 


(III) PHystcaL DISCUSSION OF RESULTS 


(1) Pressure Distribution Inside the Duct 


In linearized approximation, the pressure coefficient 
turns out to be: 


?P — Po bod 20; 


oe oa “ 
‘ (1/2)pU? U 
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Fic. 9. Pressure distribution for constant sources. 


The pressure distribution is known as soon as 2; is de- 
termined. 

(a) Constant Sources——In Case (A) of constant 
sources, one obtains, by differentiating the perturbation 
potential Eq. (5) with respect to x: 


— 
Vv, = -2ad f du(x? — B?R?)~"? (36) 

0 
Eq. (36) turns out to be essentially a complete elliptic 


integral of the first kind, A. After some transforma- 
tions, one obtains: 


6 = 48-*¢-!9 Y~ K(k) (37a) 
for the front district: 1—- Y <Z<1+ Y; and 
cy = 48-91 Yh, K(k) (37b) 
for the rear district: 1+ Y < Z. 
Here again: 
k? = (Z? — (1 — Y)?J/4Y 
k,? = 4Y/(Z? — (1 — Y)?] 
On the axis of the skeleton cylinder (Y = 0), one 
obtains from Eq. (36) for the pressure coefficient: 
Coe = 48-%9(Z? — 1)7"” (38) 


This holds for the rear region, the front region being 
restricted to a single point Z = 1, Y = 0. 
From Eq. (37a) we see that, on the front cone (k = 0), 


oa om 28-°e yoys (39) 


A discontinuity in the pressure field will develop which 
represents a conical compression shock that is weak ex- 
cept in the immediate vicinity of the vertex. A com- 
parison of the linearized solution with the solution of 
the exact differential equation is given at the end of this 
section. 

On the second Mach cone, k = k; = 1; then the ellip- 
tic integrals and c, become logarithmically infinite. 
This singularity of the pressure field is due to the fact 
that the source intensity increases at the entrance rim 


discontinuously from zero to the value A, corresponding 
to a finite lip angle, and that this discontinuity is prop- 
Although the 
prerequisites of the linearized theory are not fulfilled 


agated along the characteristic surface. 


here, it is explained in Section ITI, 2, that the region to 
be excluded from quantitative interpretation is only 
small. 

The graphs of Fig. 9, showing c,(Z, }'), plotted for 
Y = 0, '/s, and 1, give insight into the variation of 
pressure coefficient inside the duct. The study of this 
case is an instructive intermediate step before consider- 
ing the c, field due to linear sources and that due to the 
source sink combination described before. One may 
also note that, in spite of the singularities of the velocity 
and pressure field, the stream function y remains regu- 
lar, so that the flux balance of the flow pattern is pre- 
served. 

(b) Linear Sources.—In the case of sources with 
linearly increasing intensity the expression for the 
pressure coefficient is easily obtained by use of the re- 
currence theorem proved in Appendix II. Since, by 
virtue of this theorem: 


, , Z 
lin “a o const. o rzp_ CONSt./ ry 
Ur = (Bd) . dZ vy (Z’) 
Co o 1—¥} 


one obtains from Eqs. (37a) and (37b): 
°Z 
: oe ; te 
— = 48-7 lg’ J} fs K(k) dZ 
J i-—} 


for the front region, and 
C lin = nas (1 4 y. 1’) s, 


Dp ae 
4B-*9 ty f ki K(k) dZ’ (A40b) 
1+¥ 
for the rear region, k and k, having the same meaning 
as previously. 

It is clear that, c,'"" being essentially an integral over 
a function that is finite everywhere except at Z = 
1 + Y, where it admits a logarithmic singularity, will 
remain finite for all finite Z-values. The difficulty 
arising in the case of constant sources is, therefore, re- 
moved in the present case. 

It is not possible to reduce integrals (40a) or (40b) 
into a combination of tabulated functions; therefore, 
approximate methods must be used to evaluate c,'" at 
an arbitrary test point. A series expansion for c, is 
available, but, because it converges slowly, graphical 
integration has been preferred. 

Still, since the integrand becomes infinite on the sec- 
ond Mach cone (Z = 1 + Y), such procedure is ap- 
plicable only outside the interval (1+ Y-—-ypw<Z< 
1 + Y + wp), uw being small compared with (1 + Y). 
To obtain the proper junction, at least two additional 
values of c, must be evaluated analytically, say c,(1 + 
Y, Y) andc,(1 + Y + u, Y). The method used for 
evaluating these two values of c, is derived in Appendix 
ITT. 
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Proceeding in the outlined manner, graphs of c,""(Z, 
Y) were constructed for Y = 0, !/2, and 1, as shown in 
Fig. 10. It may be noted that each curve (except that 
for Y = 0) rises gradually with increasing Z, being’ at 
first concave and then convex with regard to the Z- 
axis. On the second Mach cone the curves possess an 
inflexion point, admitting a vertical tangent—i.e., a 
sudden pressure increment, as in a compression wave. 
For large Z, all lines are seen to become indistinguish- 
able, and for Z — o they tend logarithmically to in- 
finity. 

Again, the case of sources with linearly increasing 
intensity does not lead to a physically admissible flow 
pattern. Although c, is finite at all finite points of the 
Z-Y-plane, it becomes infinite for large values of Z, 
which means that the area of the exit section of the 
stream tubes tends to zero, resulting in an infinite 
drag. 

(c) Source Sink Combination.—We are now in a 
position to determine under what conditions an ad- 
missible flow pattern may be obtained. Let the source 
distribution f(¢) on the skeleton cylinder be an arbi- 
trary function of ¢, provided that f(¢) is continuous 
everwhere and that f(0) = 0. Such a function may be 
approximated by a superposition of functions of the 
type: 
for F > f, 
for ¢ < &, 


Id) _ A,'(¢ - Cu), 
flO) = 9, 


Let Cy,(Z) 


where the ¢, form an increasing sequence. 
be the corresponding pressure field. If 


\ 
A) = DV fuls) 


u=1 


then 


\ 
E(Z) = >. Gn(Z) 
n=1 
where ¢, is the pressure distribution corresponding to the 
source distribution f(¢). Since the c,, are finite every- 
where except for Z — , ¢, is likewise finite. For Z— 
©, one obtains, by using the series (32) for (1/Y) 
(0y""/OV), which according to Eq. (12b) is proportional 


-, lin, 
tov, ° 


N 
c= 48-2 } a, log 2(Z — ¢,) — 
n=1 
N 
(4 Drops + (2) 
4 : (Z a f.)* 32 
N 
o , 
1 1y2 y3 ) ——*—. 4 .,, 4 
(1 + + ) (Z — f)4 - (41) 


n=1 
with 


on’ = A,’ Br/U 








1 ee 3 


Fic. 10. Pressure distribution for linearly increasing sources. 


A more involved argument shows that ¢, finally takes 
the form: 





where fs(¥) are polynomials in Y of order 2s, which re- 
main bounded for Y < lass— ~. 


From Eq. (42) it is evident that, for Z — ©, c, will 
N N 


remain finite only if > o,/ = 0—-ie., > A,’ = 0. 


n=1 n=1 
If this is the case, then 
lin ¢(Z, Y) = 0 
Z— © 

The meaning of the above is that, in order to eliminate 
all the singularities that may eventually arise in the 
flow pattern, the source intensity function f(¢) has 
to be selected so as to approach a constant value for 

large ¢-values and to vanish for ¢ = 0. 
As an instructive application, one may consider again 
the superposition of fwo linear source distributions— 


namely, 


N= Pa A,’ = —A,! = Pi 
£& = 0; & = B 


Putting B/8d = 0.25 and 2.0, respectively, graphs of 
cy(Z) for Y = 0, 1/2, and 1 are shown in Fig. 11. One 
can notice a relation between the form of the wall con- 
tour and the pressure distribution along the Z-direc- 
tion. Beginning at the rim, there starts naturally a 
slow compression, indicated by a rising c,; in the region 
where the duct wall becomes convex toward the axis 
(front district of sinks), the c, rise is almost annulled. 
At the second Mach cone of the sources (shortly be- 
hind the throat), a sudden pressure increment occurs, 
indicating the formation of a compression wave due to 
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Fic. 11. Pressure distribution for source sink combinations. 
the combined effect of all parts of the rim. At the 
second Mach cone of the sinks (near the inflexion point 
of the wall contour), a sudden pressure drop occurs, 
indicating the formation of an expansion wave. 

(d) The Drag._-Application of the momentum theo- 
rem to the fluid in the interior of the duct yields, for the 
drag: 


D = Ao(po + poVo?) — Ai(pi + piV1") 


where V stands for the total velocity and the subscripts 
0 and | refer to the entrance and exit sections of the 
duct. It had been shown previously that for Case (A) 
of constant sources and in Case (C) of the source sink 
combination Ay = A, Vo = Vi, and pp = pp. 

The discussion of the last paragraph has shown that 
for Cases (A) and (C): c¢,(~7) = 0O—i.e., po = px = 0. 
Therefore, D = 0. 

As mentioned before, in the Case (B) of linear sources, 
D= o, 


(2) Constant Source Distribution as Limiting Case of 

Source Sink Combination 

As mentioned in the Introduction, the velocity and 
pressure field, produced by a distribution of constant 
sources (A), may be considered as limiting case of the 
more general flow pattern produced by the source sink 
combination (C). Qualitatively, this might be anti- 
cipated by comparing the pressure curves in Figs. 9 
and1l. Upon reducing the distance B, between sources 
and sinks, gradually to zero, the pressure increments 
indicating the compression and subsequent expansion 
waves, which are at distance B in Fig. 11, approach 
each other, until, in the limit B — 0, described by 
Fig. 9, the expansion wave follows the compression 
wave immediately without interval. During this limit- 
ing process the pressure values between the two waves 
increase gradually beyond any bound; this indicates, 
on the other hand, that these c, values obtained during 
the limiting process are outside of the linearized theory. 
However, if the limiting process is not carried to the 
end, the c, values remain permissible within the linear- 


ized theory and indicate approximately the true form 
of the pressure curves. 

Consider once more the source distribution (C), rep- 
resented in Fig. 12 and resolved into a combination of 
sources and sinks that are a distance B apart from each 
other. While reducing B but keeping (B/8d)A’ = 
A = constant, A’ must increase indefinitely; thus, for 
sufficiently small B values, the intensities +A’ of the 
linear sources and sinks (i.e., the slopes of the two 
parallel lines) will become so large that the prerequisites 
are no longer fulfilled for applying the linearized theory 
to test points P lying inside stripes of width B, adjacent 
to the first and second Mach cones. 

It is instructive to verify quantitatively that the 
formulas derived for the source sink combination 
[Case (C)] become, in the limit B — 0, identical with the 
formulas derived for constant sources [Case (A) }. 

Symbolically, the velocities of the source sink com- 
bination may be written in the form: 


i pilin (x _ B)] 


. li \ . lin; 
lim A,v'"(x) = lim [v'"(x) 
B-—0 B—O0O 


and after expanding in the vicinity of B = 0: 


; dvii2 
lim A, v(x) = lim | 3 -) | 
B—0 B—0 Ox 


Now (B/8d) = A/A’, so that, by use of the recurrence 
theorem, 


lin 
const 


lim Ag,o'"(x) = (8dA/A’) — =; (x) 
B—>0 Ox 

One may also carry out this check for front and rear 
district explicitly and follow up, in detail, how in the 
limit B — O the singularities at the first and second 
Mach cones gradually build up. (This is shown in 
Appendix IV.) But the true behavior of the pressure 
field is better approximated by considering the c, curve 
obtained with finite B. 

It is interesting to note that the occurrence of a com- 
pression wave followed closely by an expansion wave, as 
well as the singular behavior of the pressure field, if the 
two waves follow each other directly, also occurs in 
supersonic flow about slender bodies of revolution. 
If, following von Karman and Moore, we put linear 
sources and, a distance B behind, linear sinks on the 








WwW, JA 
B-O 


Case (A) as limit of Case C. 
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axis of an ogive, the pressure field shows also in this case 
a pressure rise, followed by a pressure drop, a distance 
B thereafter. Again, if distance B is gradually reduced 
to zero, the pressure increments indicating the com- 
pression and expansion, approach each other and grow 
beyond any bound. In this case, the conical surface, 
originating at the tip of the ogive and representing an 
attached shock wave, appears in the linearized theory 
as singularity of the flow field, at which c, becomes in- 
finite as /V1 — (8r/x)*. Again, the true behavior 
of the pressure field is better approximated by con- 
sidering the c, curve obtained with finite B. 


(3) Computation of the Form of Entrance Shock, in Second 

Approximation 

In the preceeding section the singularities of the 
flow field, produced by constant sources have been 
shown to build up gradually during a limiting process 
that allows to interpret actual conditions in a qualita- 
tive manner. In the present consideration, an at- 
tempt will be made to give further elucidation by com- 
puting the form of the entrance shock in second ap- 
proximation and by comparing these results with 
Ferri’s® analysis, carried out by means of the method 
of characteristics. Among other things, Ferri has 
shown that, for small lip angles 6, there is a front shock 
and a reflected shock, each of which has an almost 
conical shape,* as expected from the linearized theory. 
According to Fig. 8 of Ferri’s paper, a lip angle of 5° 
produces shock fronts, the curvature of which is slight 
except in the immediate vicinity of the vertex of the 
double cone, this vertex region covering approximately 
1 per cent of the total area of the shock. One can also 
read off from Fig. 8 that the velocity increment across 
the front shock follows approximately a Y~'/-law, if 
one does not approach the axis closely. 

From the linearized theory, a first idea of the curva- 
ture of the shock front may be obtained by inserting the 
velocity field obtained heretofore into the equations of 
the shock polar.’ These equations express the slope 
of the shock front \ in terms of the velocity components 
before and behind the shock front: 
uy = (1 — p?)m, cos? X + (cx?/m,) ] 
= —[(1 be”), Cos? uy, + (43) 


Cx?/u,] cot X = (uy — ue)cot y) 


Vo 


here 


“uy = U 
Cx? = (U2/M?)(B%u? + 1) 
w (y -—1)/(¥+ DD 


For u2 and v, the velocity components behind the shock, 
we substitute the values, at the Mach lines, as obtained 


from our linearized theory: 

* Interferometric pictures of a supersonic jet, emerging from a 
circular nozzle, taken by R. Ladenburg,® show clearly com- 
pression and expansion waves forming a pointed double cone. 
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Fic. 13. Shock form, in second approximation, 
TA 6 
uy =v, = U — ,= u(t — =) 
BY } BY Y 
TA 6U 
ve = 0, = = 


Vy vY 
here, 6 is again the lip angle of the duct, which, ac- 
cording to Eq. (24), isa/8. 

The solution of either one of Eqs. (43), yields, up to 
the first power of 6 (6 considered small) : 


tan A = 


6M ) dY 
dZ 


3( 283(1 — w)V¥, 


which, when integrated, yields as equation for the ap- 


proximate shock line: 


5M‘ 
Z=(1— ¥Y)- di- Y) (45) 
287(1 — mu?) © ¥ 
in qualitative agreement with the exact theory. (See 


Fig. 13.) It is evident that this result is valid only in a 
region not too far from the nozzle wall, since in the solu- 
tion procedure use was made of the smallness of 5—i.e., 
the smallness of #2 and v. In the immediate neighbor- 
hood of the wall, the expression (44) for A can even be 
shown to be a solution, if one uses the exact differential 
equation. 

Proceeding further along tiis line, higher approxi- 
mations of the velocity and pressure field, near the 
shock, can be computed. With these corrections, it 
can be shown that the singular values of v and c,, at the 


axis, disappear. 


(4) Comparison with Measurements 


It is, of course, desirable to compare the results of the 
theory with available measurements. Obviously, such 
comparison can be carried out more successfully after 
the ‘‘inverse problem” has been solved so that the pres- 
sure field of an experimentally given nozzle shape can 
be computed. This will be done in a separate paper 


(see dissertation of Ch. Mack). 
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APPENDIX 1: PROOF SHOWING THAT EQ. (2) FULFILLS 
THE AXIAL SYMMETRICAL POTENTIAL EQUATION FOR 
SUPERSONIC FLow 


It is known that the Laplace equation (in cylindrical 
coordinates) is fulfilled by the velocity potential of an 


incompressible source : 
p~ = c/V (x — §)? + R? (Al) 


where c is the source strength, (v — &) is the axial dis- 
tance, and R is the radial distance between P and the 
source element, as before. Suppose that incompres- 
sible sources of constant intensity c are spread over a 
cylindrical surface of axial length Z, which is sub- 
divided into two parts: a cylinder covered by sources 
0 < — < x, and another one covered by sources x < 
~ <L. The contributions of the two cylinders to the 
potential at P(x, r) shall be considered separately. The 
solution, representing the contribution of the first 
cylinder, assumes the form: 


"wy “x f(£) 
op = 2d f dw / dé 7 A2 
0 J 0 Vix — &)°? +R 


This expression will certainly fulfill the axial svimmet 


r 


rical form of the Laplace equation: 
Ob ob | Of 
os Or r or 


0 AS 


From the point of view of the theory of complex func 
tions, x and £ may assume complex values and still 
fulfill differential Eq. (A8). Therefore, if we substitute 

x— & = (1/#8)(x’ — & (A4) 
where x’ and ’ are real variables, then we get a solu 


tion of the form: 


mts ol f(é) ' 
@ = —2d dw dt’ = (Ad 
0 J/ 0 Vv (Xx * ‘2 B°R 


obeying the differential equation, 


fee 


ob OO 1 


Bg? — . - () (A6) 
Ox'* oF* r Or 


which turns out to be the equation determining the 
supersonic velocity potential. 

It is easily seen that Eq. (A5) becomes imaginary for 
t’ > x’ — BR. Therefore, the value of ® will be com- 
plex. Since the real, as well as the imaginary, part of 
® has to satisfy the differential equation separately, one 
may conclude that Eq. (A6) is fulfilled also by: 


We *x—B8R ‘fc ( g ) 
pm = “ 2d dw / dé J 
0 J 0 V(x — £)? — BR’ 


which is the desired, real value of the velocity potential. 
The contributions of the second cylinder covered by 


sources x < & < L to # is zero, since sources down- 


(A7) 
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stream of P cannot possibly influence © according to 
the laws of supersonic flow. 


APPENDIX II: A RECURRENCE THEOREM, OF USE IN 
SUPERSONIC SOURCE DISTRIBUTIONS ON CYLINDRICAL 
SURFACES 
Let ®; be the potential at the point (x, 7) due to a 
distribution of sources or sinks of local intensity f(¢ 

£ > Oona skeleton cylinder, viz. : 


Wx "x—BR f(£) 
®, oa f aw f dé F = 
J 0 a © V (x wie £)? —_ B?R? 


where R? = d? + 7° 
(x, 7) is located in the front district, or wx = 7, if (x, 7) is 
#, is a solution of the par 


AS 


2dr cos w, and R(wx) = x/8, if 


located in the rear district. 
tial differential equation: 
0" 


O° 1 Of 
+ O° _ 
Ox" 


or y or 


0 (AY 


If we put 


then 
P, = Id / dwF(x, R) 
J 0 
Integrating Eq. (A10) by parts, we obtain: 
F(x, R) = 
/ a = 
) f(¢) log l(a — 4) +V (x t)? —B?R?| > 
t=0 
**—BR 
eae om 
/ dé f’(&) log I(x £) + V (x — £)?— BR I 
or 
Fix, R f(O) log (xv + Vx? B?R?) — 
*x—8R 
f(x BR) log (BR) + / dé f’(£) log 
J 0 
" / ’ ae oD: , 
(vw — 4) + V(x — §&)? — BR] (All 
Now, 
Ob Ow 
(v,) L = —2MFlx, R(we)} — — 
Ox Ox 
“we OF (x, R) 
2d dw (Al2 
J 0 Ox 
However, in the front district, R(wsx) = x/8 and Fix, 
R(wx)| = 0; whereas in the rear district, w. = 7 
and Ow, ‘Ox = 0; hence, 
' OF(w, R) 
(v-)} = —2d dw Al2a 
J0 Ox 
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and after substituting F(x, R) from Eq. (All): 


0) 
poof elles 
x—BR P £) 
f — ( (A13) 
) = = Pe. 


Now, consider a second source distribution g(£) re- 
lated to the first by 


€ 
g() = a/ty fe ae’ (A14) 


where L is some characteristic length. Obviously, 
e(0) = O and g’(t) = (1/L)f(£). Let &,, be the po 
tential due to this new source distribution. Then the 
expression (A13) certainly applies to g(£) and 


1.€., 


OP;; idea 
ie = (v,) II ae 


dé V(x 


which is equal to Eq. (A8), except for the factor 1/Z. 
Thus, in general, front and rear: 


O;,/0x = (vz) 41 = (1/L)#, 


st 
2 _ BR? 


(A115) 


provided that the source distribution corresponding to 
?,, is related to the source distribution corresponding to 
®, by Eq. (A14). 

Now, according to Eq. (12b) and Eq, (A15): 


, OP, 
by = —poB? J dr’ r’ — = 
YI Po f a 
y O°? ;; 
—po8°L Ir — (A16) 
ae I Ox? 


Since, however, ®,; is a solution of the partial differ 
ential Eq. (A9): 


O°; 0°?;; l OP; 1 oO ( oe 
ee =——o——— SE Y= — = 


Ox or? r Or r or or 
= ion 1ALe 
a ( (A17) 
we obtain, by substituting Eq. (A17) into (A16): 
v1 = —polr(v,) 1 (A18) 


In order to determine L, consider the case of constant 
and linear sources: Let 
f(t) = A = Uc/Br (constant sources) 
and 
g(t) = A’(t/Bd) = 


(U/Br)o'(§/Bd) (linear sources) 


Obviously, then, 


(é) = (<) P f(t’)dé’ 
g = Bd s Jo Ss Ye K 


FLOW 
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and comparing this relation with Eq. (14), 
L = Bd(a/o’) 


Hence, finally, from Eqs. (Al5) and (AIS) we ob 


tain: 


(A19) 


— pod 3 (“,) Vy," (A20) 
0 


ve" = (pd)— 


Oo ’ 
Vironst _ — pod B ( ) ry,/in = 
o 


and with Eq. (12b): 


(o' 6) Peun 


“4 
yin = (¢'/c) / GZ." Wins (A21) 
eg 
These formulas have been used in the text. 
AppeNDIX III: A METHOD FOR EVALUATING ¢,'" AT 


THE Points (1) Z = 1+ Y, Y aAnp 
(2)Z=1+ VY+u, Y, WHere p/Y <1 


(1) On the second Mach cone, a well-converging 

series expansion has been derived. It can be shown 
. ~ li 

that, after some transformations, Eq. (40a) for c,"" as 


sumes the form: 


ol + Y, ¥) = 4o’n-"8-1(1 + Y, Y) (A22) 
here, 
ae | + y cos 0 
1+ Y, Y) = die —"— Cm 
0 1 — y cos @ 
and y = 20/Y/(1 + Y). For y = 0—ie., Y = 0— 
Eq. (A23) gives immediately, 7 = 0. Observing that 


for Y 1, y = 1 and that for Y = '/2, y = 0.94, we 
introduce a new parameter e, defined by 
e=1—-—y=(1 -— VYV)?/(I + Y) 
and rewrite Eq. (A23) as 
Ti+ Y, Y) = 244+ 14+ 2d 


with 


wr/2 \ 
I, = ‘ dé log (1 + cos #) = 
0 
- ( 2% = 
ees oe 2u : 1)? 


n=0 


1.83 — (w/2) log 2 


I, = f d@ log (1 + y cos #) = 
0 | 


r log (1 + Y) 


*n/2 0s 0 
. / d@ log | ~ ( — )|- 
0 1 + cos 0 
4 n «/2 " m 
° ¢ 
- oS fae & (-y" 
“f m=0 (1 + cos 6)™ ) 


n=1 








(A24) 





where C,," = binomial coefficients = n!/m!(n — m)!. 
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The integrals involved may be computed by means 


of a recurrence relation giving 


; adi dé | m—1 . 
Sin = ee oa Tus —1 
0 (1 + cos @) 2m — 1 2m — 1 


Hence, 


Iz = —0.57le — 0.119e? — 0.03568 — ... 


and, finally, forc,(1 + Y, Y): 


1+ Y 
cpX(1 + Y, Y) = 4r—'8 *o'[x log ( : )+ 


3.68 —1.142e€ — 0.238e? — 0.070e? 
In particular: 
Y=le=0 Cyp/4B—7-2x— 10’ = 3.68 
Y = !/,,€ = 0.06 = 2.44 


= 0.00 


(2) To evaluate c,(1 + Y + u, V), closely behind the 
second Mach cone, set in Eq. (40a) and (40b) for 
K(k;): 

K(k) ~ log (4/k1’) 

change the variable of integration from Z ton = Z — 
(1 + Y), expand the integrand into a series of the type 
> (a, + Bn log uu", and integrate. This yields, ap- 
n 
proximately, 
Gti Y +2, ¥) = 61 + J, ¥) + 

pte to" Sa(1 $1 = )+ 

2p *ax~*o * <p og = [Oe pF + 

' 1+) 


1+ Y jet 2Y + 3¥? 
FPL 20 + Y)? 
32Y 


I+ yY 


+ log “| + .. , 
This expression is valid for /Y <_ 1 and yields, for 
u = 0.05, Y = 1—ie., Z = 2.05 
Cp = 4x 18—o'(3.68 + 0.228 
for» = 0.05, Y = '/, 
Cp = 40-18-20 (2.77 + 0.229) 
For Y = 0, when the above formula becomes invalid, a 
closed formula may be easily obtained; using Eq. (38) 
one gets: 


C, = 48-*o’ 


48-*o' log (Z + VZ?—1) (A25) 


After the values of c, at Z = 1+ Y + uw have been 
determined analytically, the remainder of the c, 
curves can again be obtained by graphical integration. 
Since, however, such a procedure leads to a cumulative 
error and since the curves c, (Z, Y = constant) approach 


lin 


each other asymptotically for large values of Z, it is 
necessary to check their respective positions. Such 
check is readily obtained, since a series expansion for 
cy after Z~ is available. 

Using Eq. (52), we obtain: 


y,/in 2 Oy 
7 . oa ett 
13-20’ tox 2% ant y?) + sone 7 
(1+4¥2+ V*)+...] (A26 


Hence, for large values of Z, the c,-curves become 


asymptotic with Eq. (A25) and 


CAZ, UV) — 6fZ, F") 


3 <0 
y — } 
for Z large enough. By means of Eq. (A26), the fol- 
lowing values were found for Z = 10: 
y=0 cy" /4a—'8-2a’ = 9.4035 
= 0.5 = 9, 4015(9.433 
1.0 = 9.3955(9.423 


For comparison, the figures obtained by graphical inte 


gration are placed in parentheses. 


APPENDIX IV: PROOF SHOWING THAT THE VELOCITIES 

AT FIRST AND SECOND Macu CONE, OBTAINED IN 

Case C, Become, In Limir B — 0, IDENTICAL WITH 
VELOCITIES OBTAINED IN CASE A 


Using the front representations [Eq. (30)] of sources 
and sinks {district (3)|, one obtains: 
Azv,""(x) = 2,/"(x) — 9,/"\x — B) = 

—28-'Ue'UY'{s[1 — a(Y)s — @(Y)s?—... - 


s’T1 — a(Y)s’ — @(Y)s’? —..... 3} 


= — 28-1[e'V “Xs — gs’) x 
fl — o(YV)(s + s’) — cof Y)(s? + 2Qss’ + st) — 


Now 


; B A 
§ S = As = == ~ 
B(2r) A’(2) ) 


so that 


Azo, '™(s) = As Az,/"(s)/As = 
—AgwY™' 1 —a(Y)(s+s’) +... 


and 
. 1 y—l1/s ce st 
lim A,zv,'"(s) ~ —A 2-2) 
BO 


As long as B is not yet zero, district (3), in which *he 
above representation holds, does not contain the point 
Y=0,Z2=1 )’, and no singularity in the flow field 
arises. 

In the vicinity of the second Mach cone, using the 
rear district representations (32) [district (6)], one may 
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it is check on how, in the limit B — 0, the logarithmic singu- ; . o ra 
such larity of the velocity field gradually builds up: im As — log [Z? — (1 + Y)?] 
| lor . * 4 
Awe) A / ° dp | | As long as B is not yet zero, one has to take the dif- 
y pb p LV (a — p)(p — B) ference quotient, As(Av''" As), instead of the deriva- 
(Wem 5 wa? on p) tive (A A’)(dv'"'"/ds); in this case, the above represen- 
a ' : : ; tation, which is only valid for test points lying in the 
where c = Z’andc’ = (Z — Zp)’. After expanding, rear districts of both sources and sinks, does not include 
a A' fd “a =| l | any position of P, for which Z = 1 + Y. In the limit 
Ani Vy (“) / » LV = 06 ~ 6 Xx B - > 0, ¢ is allowed to become equal to a, and the eXx- 
7. q pression in the denominator of the integrand indicates 
= Vc—p+.... that a logarithmic singularity will be formed. But the 
dc true behavior of the pressure field is better approximated 
me - (72 ) Z(bd)— [ dp x by considering the c, curve obtained with finite B. 
Jb p 
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Two-Dimensional Airfoils in Hypersonic 
Flows 


RICHARD D. LINNELL? 
Massachusetts Institute of Technology 


SUMMARY 

Equations for the force and moment coefficients of a sharp 
edged two-dimensional airfoil aredeveloped on the basis of Tsien’s 
hypersonic similarity law. The coefficients are calculated for a 
flat plate and for three airfoils with thickness. The optimum air 
foil is a wedge whose lower surface is flat, with the maximum 
thickness at mid-chord. The results are valid when the product 
of the Mach Number and the thickness ratio is of the order of 1 


INTRODUCTION 


oe FLOWS are flow fields where the fluid 
velocity is much larger than the velocity of sound. 
Hypersonic flows have been treated qualitatively, 
usually on the basis of the limiting case of an infinite 
value of the free-stream Mach Number, 1/°, which is 
defined as the ratio of the undisturbed fluid velocity to 
the velocity of propagation of small disturbances, the 
velocity of sound. Von Karman! has considered this 
case when y, the ratio of specific heats, is 1, for which the 
relations are particularly simple. He points out that 
in many ways the dynamics of hypersonic flows is 
similar to Newton's corpuscular theory of aerody- 
namics. The pressure acting on an inclined surface is 
thus greater than the free-stream pressure by a quantity 
that is approximately proportional to the square of the 
angle of inclination, instead of the usual linear law of 
conventional flows. Sanger’ has used this concept to 
design the best wing and body shapes for flight at an 
infinite value of M°. The qualitative approach, how- 
ever, is not satisfactory for hypersonic flows whose free- 
stream Mach Number is finite. 

Tsien*® has developed the similarity laws for hyper- 
sonic flows. An affine transformation which expands 
the flow field laterally, reduces the equations of the 
flows to a single nondimensional equation. If a series 
of bodies having the same thickness distribution but 
different small thickness ratios, 5, are put into different 
flows such that the parameter K = 11% remains con 
stant, and at angles of attack, a, such that the ratio a/6 
remains constant, then the flows are similar in that they 
are governed by the same transformed velocity poten- 
tial. A must be of the order of 1. 

The hypersonic similarity laws state that the force 
and moment coefficients for two-dimensional flows can 
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* The material in this paper is taken from a Master’s thesis at 
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be expressed in the forms 
C, = @L(K) 


Ca = &D(K) 
Cn = 6°M(K) 


The purpose of this paper is to determine the func 
tions L, D, and M/ for two-dimensional airfoils. 
THE PRESSURE COEFFICIENT 
According to the hypersonic similarity law, the pres 
sure coefficient on the surface of an airfoil in a two 
dimensional flow can be written in the form 


C, = &P(K, £) 


Here 6 is the thickness ratio of the body, the ratio of the 


maximum thickness of the body to its length; A is a 
parameter defined by K = 4%; and é, a nondimen 
sional coordinate along the chord, is defined bv the 
affine transformation 


x bé 


y= bin 


In the x, y axes the midpoint of the 
Its chord, of length 2), 


of the fluid field. 
body length is at the origin. 
is on the x-axis. 

The pressure at a point on the surface of the airfoil 
will depend on the free-stream pressure p°, on the pres- 
sure change due to an oblique shock attached to the 
leading edge to a pressure p,, and on the pressure change 
through isentropic expansion to the point on the airfoil 
and a pressure p. Thus, two relations (p,/p°)(11%O,) 
and (p/p,)(.\/°6,, 17°0,), must be determined. 0, is the 
deflection angle produced by the shock, and @, is the 
angle through which the flow expands. 

4\7°6, enters in the expansion pressure ratio through 
the Mach Number after the shock, W,, which will be 
the initial Mach Number of the expansion. , is de- 
termined by the conditions of the shock, so M,(1/°9,) 
must be found. 

6, is a function of the angle of attack and of the slope 
of the body at the leading edge. It does not vary with 
£. 6,, however, is a function of the angle of attack and 
of the thickness distribution along the chord. It is then 
a function of §. After these relations, in addition to the 
pressure ratios, have been found, the pressure coefficient 
can be written in terms of K and £. 
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TWO-DIMENSIONAI 


[THE SHOCK RELATIONS 


lo determine the shock relations consider an oblique 
attached shock at the nose of a wedge of half angle @,. 
a is the angle of the shock relative to the free-stream 
flow direction. /p,, the pressure immediately behind the 
shock, will also be the pressure on the surface of the 
wedge. The flow will be deflected through the angle @, 
by the shock. 

By definition, the pressure coefficient is 


car 23 (? 4 
1/2)p°U y p? U=\p* > 


Thus, 
p,/p” (1/2)yM°C,, + 1 (2) 
The thickness ratio of the wedge is @,. The flow 


parameter is then A 17°,, and from the similarity 


law 


"oe 6,.°P.(Ks) (3) 


D,/ p' (} o)yK,;?P, + ] (4) 


The function P, must now be determined. p,/p® is 
given by the relations for oblique shock as‘ 





Ds 2 ee y= iI - 
= M° sin? a — (5) 
- vt} ++1 
Substitution of Eq. (5) in Eq. (1) gives 
Cy. = [4/(y + 1)][sin? a — (1/M°’)] (6) 
or 
sin? a — (1/M®’) = [(y + 1)/4]CGs (7) 


In order to eliminate a, a relation between a and 6, is 


needed. For this, consider‘ 
1 ie y + 1 sin asin 6, , 
; = sin’ a - ; (8) 
if 2 cos (a — §,) 
which, in consideration of Eq. (7), becomes 
sin a sin 6,/cos (a — 0,) = ('/2)C,; (9) 


Finally, by expanding the cosine term, Eq. (9) becomes 
('/,)C,,2 cos? a = [1 — ('/2)C,,]* sin? a tan? 6, (10) 


Eq. (10) is still general. It is now necessary to make 


the approximations appropriate for the conditions of 
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hypersonic flows. For small values of @,—that is, for 
thin airfoils at small angles of attack—C,, will be small. 
For the condition that .1/° is large, in conjunction with 
small 4,, a will also be small. Therefore, the appropri- 
ate approximations are 


[1 — €/2)Cp? 1] 
cos*° a => ] (11) 
tan? 0, = 0.24 


With these approximations and Eq. (7), Eq. (10) be- 
comes 


: ee me ; I ‘ 
iC,, = 0,7 sin? a 6, , Cy. + ye: (12 


The substitution of C,, = @,°P, and K, = 11%, reduces 
Eq. (12) to 


P,? — (y + 1)P, — 4/K2 = 0 (13) 
And, finally, 
P, = [(y + 1)/2]) + V[(vy + 1)/2]? + (4/Ky?) (14) 


the positive sign of the radical being required since P, 
must be positive. 

The value of C,, when /° becomes infinite is needed 
in order to make a comparison with Sanger’s results. 


Coe > (7 + 1)8,, as M°— (15) 


In a similar manner the equation for .1/,, the Mach 
Number after the shock, is found to be 
“ef ge 
o Vi 7 2p ( x 
M ((y + 1)/4]K2P, + 1 
fl(y + 1)/4)K PP, + ik 


) (y/2)K2P, + 1 (16) 


A?P, takes on all values from zero to infinity, so the 
second term of this expression varies from 1 down to 
V6/7 = 0.92, for y = 14. 
neglected but will be included in all calculations here. 


It could probably be 


Also, by using Eq. (15), 


1f2 1 we Z 
M, 2 —{- ; as M°— » (17) 
8 Y : Je 


if K,2P, = M°°C,, is much greater than 1—that is, if 0, is 
finite. For y = 1.4, the square root in Eq. (17) is 
greater than 1; and since @, is small, 1/0, is large. 
Thus, M, is much greater than 1 in this limiting case 
This condition will be used in making the approxima- 
tions for the expansion relations that must now be de- 


termined. 


THE EXPANSION RELATIONS 


The pressure ratio for the expansion will be determined by consideration of expansion around a corner, which is 
known as Prandtl-Meyer flow. The relation between the initial Mach Number, /,, the final Mach Number, M, 
and the angle of expansion, @,, in a Prandtl-Mever expansion, is given by 
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rare 2 ee 4 
iGugit* (tan ‘ei? VM? — 1 — tan-'4/7—._ Vm, — i) s 
y¥—1 ¥+ 1 ¥+ 1 


(tan 1WM? — 1 — tan iVM,? ~ i) (18) 


For hypersonic conditions /, is much greater than 1, as shown by Eq. (17), and .V/ is greater than J/,.  There- 


fore, 


VM2—-1—=, 


VM,’ -—1—>M, (19) 


If, then, the are tangents are expanded in infinite series, Eq. (18) becomes 


1 1 1\' + 1\" 
PEs SCE)" 
y-1L2 M\y-1 3M \y - 1 


ra l (2 4 \ 
2 M,\y-1 


(7H) + ...|-(-t+ rl I 
3M,? \y —1 pes 2 M' 3M 2'M, 3M3 °° 


ae) “3 LG) ~ 
0, = ~- ao “7 _ 
y-1\"M, M 3L\7 — 1 Mi M 


) + (20) 


The relation between the pressure and Mach Number for an isentropic process such as the expansion is given by 


p  fl(vy — 1)/2)M2 + 1)" ” M\oUo~ © ie 
~ . ~ 21) 
>, Uly —1/2]M2 +15 M ' 
or 
M = M,(p/p,)~% ~ °/ (22 


Substitution of Eq. (22) in Eq. (20) gives 


9 1 p (y 1)/2y l 
ah-OY 
y—- 1M, p. 3M,3 


In order to obtain a relation between p/p, and A/,,, 
only the first term of this series can be retained. Then, 


M0. = [2/(y — 1)][1 — (p/p) 7 P77) (24) 


or 
( - 1) 97 
l (25) 


b/p, = {1 — [(y — 1)/2]M,0.}* 
If K, = M8, is defined, Eq. (25) becomes 


b/p, = (1 — [Cv — 1)/2](M/M%)K,}?"° ~ (26) 

Eq. (26) determines the value of p/p, for [(y — 1) 2] 
(M,/M°)K, < 1. For [(y — 1)/2](M,/M°)K, > 1, 
b/p;is zero. In the latter region the equation will give 
imaginary values for most values of y and positive real 
values for those values of y which make the exponent 
an even integer. 

The physical reason for this behavior is that the pres- 
sure goes to zero for a finite deflection. If the surface 
of the body deflects more than this limiting value, the 
flow leaves the surface and the pressure remains zero, 
When a given deflection is required at a point on the air- 
foil, then as 17° increases, it reaches a value at which the 
maximum physically possible deflection is required. 
The flow has expanded to an infinite Mach Number and 
to a zero pressure. This occurs when [(y — 1)/2] 
(M,/M°)K, = 1. For greater M°, the pressure at the 


given point is still zero, but the flow has separated from 
the surface farther upstream. 

By consideration of Eqs. (17) and (25), it is found 
that 


p y —1\'%0, PY ~ P 
—> —_ as 
Ds 2¥ A. ; 


M°— ~,0, #0 (27) 


The pressure coefficient can now be written in terms 


of K, and K,. 
Cy a = (2 — i) = - (2 Ps — i)- 
yM° \ p® yM°® \b, p° 


9 = 2y/(y — 1) 9 
(02. + —— (: me = K.) ~ —,, 
yu 2 MM yM° 


(28) 
For the case of no shock, @, is zero and M,/M°® is 
unity. For no expansion after shock, C, = 6,?P, 


aes 
THE DEFLECTION ANGLES 


Eq. (28) is not yet a function of K and &, so now the 
deflection angles must be determined from the boundary 
condition that the flow is tangent to the airfoil surface 
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at the surface, unless separation occurs by reason of the For values of the angle of attack less than the slope of 
flow having expanded to a zero pressure. the upper surface at the leading edge, shock will occur 
on the upper surface; for greater angles of attack, no 
The values of the deflection angles are then deter- shock will occur on the upper surface, so the angle of 
mined by the thickness distribution and angle of attack shock deflection will be zero. Thus, the deflection 
of the airfoil. The profile of the airfoil is given in x,y angles must be evaluated in the two regions separately 
axes, with the chord on the x-axis and the midpoint of for the upper surface. 
the chord at the origin. The leading edge is at — 0b, and For the purpose of writing out the deflection angles, 
the trailing edge is at b; the chord is 2b. The profile the following notation is convenient: f(x) = dy/dx = 
must be slender and pointed, as for linearized theory, slope of the profile; subscript u indicates the upper 
since the deflection angles have been considered to be surface; subscript / indicates the lower surface; sub- 


small. The x-axis is taken as inclined to the direction script —0 indicates the leading edge value; and f, = f; 
of the undisturbed flow at the angle of attack a. . for symmetrical profiles. Then, 
0 < a < f b, 65, Wie Le ~~ & G6, i f Rieu Fu(x) 
ie te 0, » = O 6.4 = a — f(x) 
all a O,, t= fu, l a Oe. ; = EZ Ae tages fi(x) (29) 
For the transformation used in the hypersonic similarity law, with dy/dt = h(£), and if a = 60, where o is the 


transformed angle of attack, the deflection angles become 


De 5 Bg: 6, » = O(h-1 y — @) 6. 4 = S5[h-1 , — h,(€)] 
orks. 6, =0 Ou = 5[o — h,(é)] 
all o 6, = b(h-1,1 + @) O..1 = b[h-i, 1 — hi(é)) 


If, finally, 6(£) is defined by 
M°e = M%8 = KB 
where the §’s are determined by Eqs. (30), then the 1/°@ parameters have been found in terms of K and &. 
The pressure coefficient can now be written in terms of A, é, and o. 
C,/26? = [('/2)B2?Ps + (1/yK*)]{1 — [(y — 1)/2](,/M% KB, }?% — ? — (1/7K?) (31) 


For a given thickness distribution, 8, depends only on the angle of attack. 6, varies with £, however. The 
evaluation of the force and moment coefficients will be simplified by the following notation, which separates those 
terms of the pressure coefficient which are not functions of £ from those that are functions of £. 


C,/252 = AB(t) — (1/yK2) 
A = [(y¥ + 1)/4]6,? + VIG + 1)/4]8.2}2 + (8,2/K2) + (1/yK?) 


B= {1 — [(y — 1)/2](M,/M%)KB,}?” — P | 
M, y¥-1, i\~"*/y —1 1 +>+1\* = 
= (? ra + A ) ¢ — —— + : } (32) 
AM? 2 2 2y? K2A 247 
0 < 0 “4 h 1, u B, = he Se B, = h = Gee h,,(&) 
0 Pg h 1 u B, = 0 B, es h wl } 
all o Bs, =hi ito Bea = hea, — hE) 


When using Eqs. (32) to evaluate the pressure coefficient on, say, the upper surface of the profile, the subscript 
must be added to all terms depending on 6’s. When writing out the full term, the correct 8’s should be used. 
FORCE AND MOMENT COEFFICIENTS 
Now that the pressure coefficient has been determined, it can be used to evaluate the lift, drag, and moment 
coefficients. The forms 


7” 5 1 I 
( 1 = (1 2b) [ ( —C,, u + G pdx = / (- 9 ts u + 5 Om as 


/-—)d ~ 


i ; . i l - 
Ca= (1/20) / [—-C,. At = Ja) 4 Cp, ila + fi) |dx = 6 / }o(- > Comm > > Cn) + 5 MtuCorn + hiCy, (33) 


b 


_ J —!1 


-} {° . : . - - i 
Se = (0p)? _(-G ut C,, )xdx = — a |, — 5 “ F 5 Cort tdé 
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TABLE | 
Coefficients for Two-Dimensional Hypersonic Flows 

= -Adet Ad Sa oe + Alt + 4d MO = At, — Ad 
0 c af + 

eae + 1 2 (7 + 1 +)". B,” ] K = M%,o = a/é 
ie 4 Br +N “i K? + 7K? M°® = free-stream Mach Number 

at ‘ . es - 

y¥—1M,,.,,, 2wvO-DN,, 5 = thickness ratio 

' rs [1 ~ 2 Me ABE | dé a = angle of attack 


VJ! 2 
Oo @ = By. Bou mh 
og a” t_ipuw Bs, = 0 
all o Bt = hi, 


h(é) 


26 on the x-axis, which is inclined at the angle of attack a@ to the flow 


When 


the terms containing it are zero. 


Sub u or /, upper or lower surface. 
terms as on @’s therein 


The same subscript appears on 
Sub —1, value at leading edge. 


—¢ Bau = hia uw m Bult 
,u =o —h,(é) 
+o 3,1 = h_yyi — Ailé) 


= (1/5)(dy/dx)(), where (dy/dx)(x) is for the surface of a slender pointed airfoil whose leading edge is at x 


For K greater than 10 the following approximations can be made 


oo eee rn) = CS | 
3g Be 2 (in)A 24 ) a 


Region of validity of equations: \/° = 5, K order of 1 


will be considered sufliciently accurate for the thin profiles and small angles of attack under consideration. 


tive moments are taken as tending to raise the leading edge 


32), 


the coefficients become 


- " 
-A, | B,dét + A, / Bydé AS. 462 
J-1 

i oe “a cd a | 
a + a, / h, Bult + A, f h Bult — =f 
6 6° . ' ; ; yK , 


= “ 
C..1.)2/o? = A, | ('/)EB,déE A, | (EB dé A 
1 1 


The last term in the drag expression drops out, since 


“a ai a 
/ h dé = / (dn /dé)dé / dn = 0 (35 
—1 al | e l 


by definition of the x-axis. The terms I, IT, and III are 
defined by Eqs. (54). 

The combination of Eqs. (32) and (34) provides a 
method for the determination of the force and moment 
coefficients of two-dimensional airfoils in hypersonic 
This method is summarized in Table 1. 


flows. 

The physical significance of the various terms in 
Table 1, especially I, II, III, becomes clearer when it 1s 
noted that A is equal to (1/yA*)(p,/p°) and therefore 
represents the pressure ratio due to the shock at the 





ne 


= —b, with a chord 


For symmetry, h, = hy. 


1) /2](.M,/M°)KB,.(é) > 1 


8,#0 


Px SI 


By substitution of the pressure coefficient from Eqs 


= 
h,d& =o 1 + A, / h,B dé 
il . l (54 
- C 
A, f hBdt = 0 “+ AM, + All, 
. l 0 
I AITI 


leading edge. Similarly, the term B of Eq. (82) is 


equal to p/p, and thus represents the expansion pressure 
ratio at any point on the chord. Thus, I, the integral 
sum of the product of the expansion pressure ratio and 
the element of the nondimensional chord, dé, over 
which it acts, is a contribution to the lift force. When 
this product is multiplied by the local nondimensional 
inclination of the airfoil surface to its chord, h(£), a 
contribution to the drag is obtained (term IT), or if it is 
multiplied by the local coordinates, £, a moment con 
tribution results (term III). 

The 8 terms clearly represent the variation of the 
nondimensional inclination of the surface of the airfoil 


to the free stream. 
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AIRFOILS 


Ww 
~“ 


IN HYPERSONIC FLOWS 


THE FLAT PLATE 


The method of determining the lift, drag, and moment coefficients of two-dimensional airfoils in hypersonic flows, 


which is summarized in Table 1, can be applied to a flat plate as a simple illustration. 


For this case 6 cancels out, 


leaving the new parameter /°a, which will be defined as Ko. 
For a flat plate h, = h, = 0. Then, for all o—that is, for all a—Table | gives the following results. 


Bs, u = 0 Bs, H = 7 
Bes u ee Be: = 0 
(M,/M°*), = 1 (M,/M°), = not needed 
I, = 2{1 — [(y — 1)/2]Ko}?”0-” I; =2 
Il, = 0 IT, = 0 
ITI, = 0 III, = ( 

1. 472 ave 3 +1\2 l l 
z = 1/y7K? A, | + C ) | 

| \ 4 + K°e? id 7K? 

Since a = 60 and K = 1/%, then Aco = Ko, so the 


coefficients are 


When ('/2)(y + 1)Ao => 1, the second term of the lift 
coefficient becomes merely 2/y)1/°. C,/a* can be 


written as a function of Ao only. 


| ee ae | } 
| M°=0 
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Lift coefficient C; of a flat plate vs. angle of attack 
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C, and C,/a? and their derivatives are continuous at 
("/o)(y + 1)Ko = 1. 

These equations confirm the statement that for large 
values of A/° the lift is approximately proportional to 
the square of the angle of attack. But they also pro- 
vide a quantitative evaluation of the degree of validity 
of this statement. It is not strictly true until 1/° is in 


finite, when 


C, = (yy 4 l)a? 


For finite values of 1/°, the lift curve starts from the 
origin as a straight line, which rapidly becomes para 
bolic as the quadratic term becomes important. The 
point of transition from a linear relation approaches the 
origin as .\/° approaches infinity. 

The values of the lift coefficient given by this hyper 
sonic equation, with y =1.4, are shown in Fig. 1 for 
f° 3,5, 10, and ©. 
\/° increases is indicated, but the most conspicuous fea 
ture is the extremely small lift coefficients available in 
The curves for 7° 0, with 


The decreasing linear range as 


the hypersonic region. 
the incompressible fluid theory slope of 27 per radian, 
and for A/° 1.5, 2, and 3, with the linearized super- 
sonic theory slope of 4/ VM — 1, are also shown in 
Fig. | for comparison. The family of curves empha 
sizes the rapid decrease of the lift coefficient as the free- 
stream Mach Number increases. For large angles of 
attack in the hypersonic region, the linear theory under 


estimates the lift and therefore is not satisfactory. 


AIRFOILS WITH THICKNESS 


The hypersonic force and moment coefficients have 
been calculated, by use of the equations of Table 1, for 
three simple airfoils with thickness. The three airfoils 


are: 
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TABLE 2 
Equations: Airfoils with Thickness 

c cS Cm(c/2) 

—y 1,1 All oY e+ AM + Ah ne 1,111, — AilIl, 

o- ) 0 o- 

Double Wedge 
—1<:<0,h4,=1,h 10<t<1k& = —-1,h4: = -1; B, i+j«,f.2" 146 
7 l g¢>1 
M, 2y l l y-l1, V2r/(y-1 

14] 1.4 Lf - y- 07). *] Ads = af 1-75 K(o — 1) | } 


1 III, - ('/)Auallu 


Lower surface equations of same form 


—1<¢€ 0, kh, = 2,4: = 0; O< 


AI, = —('/s)Aull, 
Lil, 2A, Aill, = Ailll; = 0, for ill a 


A, III, = —(}/,)A.II, 


Lower surface as for ¢ < 1 


Single Wedge 


g l,k,y = —2,4:=0: 8.4 =2-—-06,8,t =o 
o > 2 
l y¥-1  gety—D 
tl - ,[1- 5 K(o - 2) | f 
l Y ‘ : 2y/(y¥—1 
ami[1- 9 K(o + 2) | 
l y-1 s : 27/(y¥—-1) 
wl. = 2 5 [1 - 5 Ko - | a 
2 2y/( l 
=k + 5 K(o +2) | = 
1,111, = —(! i ake 


1= 0; Bau = 


{ l 2 l \ y-1 (37-1 ¥—-1) = : (8y-1)/(y-1)) 
iP 2 =- (5 mag x) \[ — 9 Kio + » | -[1 a - K(o — » | ( 


2 l \ y¥- 1 
< Al, = ( a)j[1-2 
724 3y — 1 yK3/)] 2 

( 1 A ee ee 
S\Sy — by = 3 es I 2 


All o; A, III, = —('/s)A,II. Ail, = 2A, 


(1) A chordwise symmetrical, symmetrical 
wedge, which will be called a double wedge. 

(2) A chordwise symmetrical single wedge 
lower surface is flat. It will be called a single 

(3) A chordwise symmetrical single parabolic 
whose lower surface is flat. It will be termed a 
parabolic profile. 


(4y -—-2) 
K(o + » | 


Ail; = All: = 0 


double 

upper surface at the leading edge will be called the 
whose critical angle of attack. For the three profiles consid 
wedge. 
profile respectively. 


single formed angle of attack, o, are 1, 2, and 4. 


this reason the angle of attack equal to the slope of the 


ered here, the critical angles of attack are 6, 26, and 46, 
Thus the critical values of the trans- 





The slope of the upper surface of the profile at the 
leading edge divides the angle of attack into two re- 
gions. For angles of attack less than this value, shock 
will occur on the upper surface of the profile, with a 
consequent loss of lift. Therefore it is desirable to 
operate at angles of attack greater than this value. For 


For an airfoil at an angle of attack above the critical, 


the contribution of the upper surface is increasingly 
negligible as .J/° increases. And, since most of the lift 
is furnished by the lower surface, any expansion on the 
lower surface is detrimental. 
the hypersonic region will have a flat lower surface. 


Thus, the best profile in 
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The equations obtained by use of the relations of 
Table 1, and used to plot Figs. 2-5, are shown in Table 
2. The values of A and .V/,/./°, which are functions of 


8,, are to be determined by the equations given for them 
in Table 1. 


DISCUSSION OF RESULTS 


The results of the calculation of the hypersonic co 
efficients for the three airfoils with thickness are shown 
in Figs. 2-5. Fig. 2 shows the lift curve of the single 
wedge profile; Figs. 3 and 4 present the polars for all 
three profiles, in two different scales; and Fig. 5 shows 
the moment curves for all three profiles. 

Fig. 2, the lift curve of the single wedge, is similar to 
the lift curve of the flat plate, Fig. 1. The effect of the 
thickness ratio on the characteristics of hypersonic air- 
foils can be found from any of the families of curves. 
Since K is directly proportional to both 6 and J/°, a 
change in 6 will have the same effect as a change of .1/° 
in the same ratio, and vice versa. 

In Fig. 3 the curve for the double wedge is drawn, 
with only the points for the other two profiles shown in 
the interests of clarity. The curves for A = 1, 10, and 
© fall almost on each other. The single parabolic pro 
file never has a 6L/D ratio better than the single wedge, 
but for low values of lift the single wedge is superior. A 
comparison between the single and double wedge shows 
that the single wedge is again superior, except for low 
values of lift. The transition value of the lift increases 
as K decreases. The single wedge is superior in this 
comparison. The increasing superiority of the flat 
bottom profiles as K increases is also indicated by Fig. 
3. 

Fig. 4 shows the polars to a larger scale so that the 
maximum 6L/D can be investigated. In hypersonic 
flows the wave drag is small, so the skin friction will be 
an important part of the total drag. The maximum 
5L/D indicated by Fig. 4 may be appreciably reduced 
by the skin friction. When friction is not considered, 
the single wedge has the greatest maximum 6L/D for 


K > 0.5. When friction is considered, the minimum 
value of K for which this is true will decrease. Simul- 
taneously, the. maximum 6L/D will decrease, and the 
lift coefficient at maximum 6L/D will increase. 

Fig. 5 shows that the moment for the double wedge in 
the hypersonic region is positive and increases as K 
As K increases, the center of pressure thus 
It does not 


increases. 
moves forward at a constant lift coefficient. 
reach the quarter-chord point. The moment co 
efficients of the single wedge and the single parabolic 
As K in- 
creases, they increase algebraically, as for the double 
wedge, but these profiles approach zero as a limit 
(above the critical angle of attack). 

On the basis of a smaller (negative) moment the 


profiles are negative about the mid-chord. 


single wedge profile is superior to the parabolic profile at 
all lift coefficients. This is a result of the smaller lead 
ing edge slope of the upper surface, for equal thickness 
ratios, of the single wedge profile. This difference be 
tween the two profiles also accounts for the difference in 
A criterion for 


6L/D’s noted in connection with Fig. 3. 
chordwise symmetric profiles for hypersonic flows can 
thus be made. For chordwise symmetric profiles the 
slope of the upper surface at the leading edge should be a 
minimum. The profile that meets this condition is the 
wedge upper surface. 

The above results, for an infinite 17°, have been com 
pared with Sanger’s results, and the agreement is satis 
factory. 
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Linearized Theory of Supersonic Control 
Surfaces 


P. A. LAGERSTROM* ann MARTHA E. GRAHAM* 
California Institute of Technology 


SUMMARY 


The purpose of the present paper is to make some contributions 
to the aerodynamic theory of control surfaces in supersonic flight. 
The control surfaces discussed have trapezoidal or triangular plan 
forms and are attached to an immovable surface (‘‘wing’’) so that 
a mutual interference exists between the stationary surface and 
the deflected control surface. The linearized theory is used ex- 
clusively. The basic property computed is the pressure distribu- 
tion induced by the deflection of the control surface; from this 
pressure distribution, lift, and center of pressure are found. It 
will be seen that, ip general, only a comparatively small part of 
the wing is influenced by the deflection of the control surface, so 
that most of the aerodynamic load will be carried by the control 
surface itself. However, in the cases considered in this report, 
the interference is favorable—i.e., the total lift on the wing and 
on the control surface induced by a deflection is greater than the 
total lift obtained by the same deflection of an isolated control 
surface. This depends mainly on the fact that the mutual de- 
pendence of the upper and lower surfaces near a subsonic edge 
which causes loss of lift as compared to the two-dimensional 
value) is prevented by the presence of a wing. For real fluids 
there might be some communication between upper and lower 
surfaces, because of leakage through the ‘‘gap’’ between wing 
and deflected surface. This leakage, which is not taken into 
account in the linearized theory, could cause loss of lift. There 
are many similar phenomena that tend to reduce lift which may 
not be treated in the linearized theory and which are intimately 
connected with the viscous properties of the fluid. Thus, in 
general, while linearized theory gives a fair prediction of lift and 
wave drag of a wing, it is expected to give values that are too high 
when applied to control surfaces. Correction factors will have 
to be found from experiments; from the suspected importance of 
viscous effects, it seems advisable to pay much attention to the 
Reynolds Number in such experiments. 

The present paper ts an abbreviated version of reference 1, 
which contains detailed explanations, computations, and also 
some discussion of the applicability of the theory to real fluids 
The basic formulas for pressure distribution are taken from the 
theory of conical flow and may be found in references 2 and 3 


NOTATION 
CL = lift coefficient resulting from aileron deflection 
based on aileron area 
So... = two-dimensional lift coefficient = 45m 
Cp = pressure coefficient 
CD, sim. = two-dimensional pressure coefficient = 26m 
A.R. = aspect ratio of aileron 
= mtan g({denotes position of leading edge of aileron) 
h = tan 8/tan uw (denotes position of outboard edge of 
aileron) 


= tan y/tan mw (denotes position of inboard edge of 


aileron) 


Presented at the Aerodynamics Session, Sixteenth Annual 
Meeting, I.A.S., New York, January 26-29, 1948 
Guggenheim Aeronautical Laboratory 


d = chord of aileron 

M = free-stream Mach Number 

m = 1/V M1 = tanu 

q = '/9U?2 

t = tan r/tan u 

to. = ray on which center of pressure is located 

U = free-stream velocity 

a = angle of attack 

B = angle denoting position of outboard (‘‘free’’) edge of 
aileron 

¥ = angle denoting position of inboard edge of aileron 

5 = angle of deflection of control surface with reference 
to wing 

mm = Mach angle 

p = density 

Tr = angle denoting arbitrary position on aileron or 


adjacent wing 
¢ angle of sweepback of leading edge of aileron 


> 
il 


RECTANGULAR PLAN FORM 


dbs ELEMENTARY EXAMPLE is well adapted to 
illustrate the fundamental principles. Consider a 
rectangular wing ABCD (Fig. 1) with a rectangular 
aileron AIEF deflected so that the increase in angle of 
attack as compared with zero deflection is 6. 

It follows from the superposition principle‘ that the 
effect of deflecting the wing is obtained by superimpos- 
ing the solution for (b) in Fig. 1 on the solution for the 
wing with undeflected aileron. The solution for (b) is 
the same as that for a rectangular wing of plan form 
AIEF except in the Mach cone from E. This is exactly 
the region where the interference between the wing and 
the aileron takes place. The presence of the flat plate 
at zero angle of attack to the left of EF corresponds to 
the presence of the wing. It means that the deflected 
surface may induce no downwash on EFG but that, on 
the other hand, a pressure differential can be carried 
here by the plate. Since a = 0, this pressure differ- 
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Fic. 1. Rectangular wing with rectangular aileron. 
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ential is zero outside the Mach cone from E. Without 
knowledge of the actual pressure distribution, one may 
conclude immediately that the lift coefficient in EGH 
based on the area of the deflected part (EFH) has the 
two-dimensional value 46m. This is so, for, if a new 
wing is formed by interchanging right and left in Fig. 1b 
and the two wings are superimposed, the pressure dis- 
tribution in the Mach cone from E will be two-dimen- 
sional. 

Actually, the pressure distribution may be found from 
a well-known solution by the following procedure. To 
fix the ideas, assume 6 >0. Consider a nonlifting wing 
of plan form AIEF with a wedge profile of half angle 6. 
It is well known that the pressure in the Mach cone is 
given by the formula 


Co! Conn = (1/2) are cos (—! (1) 


/ 


Since this nonlifting wing does not induce any down- 
wash in the plane of the wing, a flat plate at zero angle 
of attack may be inserted there without changing the 
flow characteristics. Assume that this flat plate covers 
atleast EFG. Then communication between the upper 
and lower surfaces is excluded in the triangle EGH. 
Consequently, if the lower surface is kept intact and the 
upper surface changed so that the wing has zero thick- 
ness, Eq. (1) is still valid for the lower surface. But 
the wing so obtained is exactly the wing in Fig. Ib. 
Thus, conditions on the lower surface are known. For 
zero thickness the perturbation pressure is antisym- 
metric in the plane of the wing—1.e., on the upper sur- 
face it has the same absolute value but opposite sign as 
compared with corresponding points on the lower sur- 
face. Since C, is referred to Cp, .4;,,., Eq. (1) is valid for 
both upper and lower surfaces in region EGH of Fig. 1b. 
Summarizing these results for the pressure distribution 
on the lifting surface in Fig. lb gives the following 
formulas. 

In the deflected part outside the Mach cones from E 
and I: 


Co/Cuan. = 1 (2a) 

In the Mach cone from I 
Cy/Cor-dim. = (1/2) arc cos (1 + 2f) (2b) 

In the Mach cone from E 


C IC 
Dp Po-dim. ~ 


; (1/7) are cos (—12) (2c) 
By integrating Eq. (2c) it is found that the loss in lift 
(as compared to the two-dimensional value) in the 


region EFH is given by 
Loss in lift in EFH = 2g6m*d?/r (3) 


This is compensated for by an induced lift in EFG of 
exactly the same amount. 

It should be noted that the plan form of the wing is 
entirely irrelevant except as it affects the area of the 
portion of the wing inside the Mach cone from E. 
Thus, a change in chord length of the wing or in the 


sweepback of the leading edge has no effect on the lift of 
the control surface according to the linearized theory. 
For a real fluid, however, the chord length will influence 
the thickness of the boundary layer and thus, presum- 
ably, also the effectiveness of the aileron. 

Also note that Eq. (2c) predicts no pressure equali- 
zation due to leakage through the gap EF. How 
much leakage there is depends on many factors beyond 
the scope of linearized theory. 


EFFECT OF SWEEPBACK 


As an example of the effect of sweepback consider the 


wing-aileron combination in Fig. 2. 





GFH \ 


Fic. 2. 


By methods similar to those indicated in the previous 
section, one obtains the following formulas for the 
pressure and lift in EGH. 


ke l (< —f ) 
: - / - are Cos : (4 
C; rV1— a l1—at 


22-dim. 

for —1<t<l,a<1—(m/A.R.). The lift coefficient 
in EGH based on the total area EGH is one-half of the 
sweepback value. However, for g > 0, the area EFH 
is greater than half the area EGH, so that one will not 
obtain the correct value of the lift by assuming that C, 
is zero in EGF and has its sweepback value in EHF. 
Actually, one obtains 


2gim*d? | : a are cos a) 


Lift in EGF = : - 
(1 — a?) \ ft ow at 9 


where a < 1 (m/A. R.). 


TRIANGULAR PLAN FORM: FREE EDGE SUPERSONIC 


Consider the wing-aileron combination in Fig. 3. 
The edge FG is referred to as the free edge. It is 
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supersonic if 8> yw. If the other edge IF is supersonic, 
there is no interference between wing and aileron, and 
the aileron may be treated as a separate wing. Hence, 
it will be assumed that |y, <u. Then there is interfer- 
ence between the wing and the aileron in the triangle 
FJK. This interference effect is easily computed by 
the method mentioned in Section 1. One starts with a 
nonlifting wing of plan form IFG and wedge profile with 
By insertion of a flat plate at zero angle 


half angle 6. 
the upper and lower surfaces are made 


of attack in IJF, 
independent, and one may change over to zero thickness 
without disturbing the flow below the wing (in the case 


§>0). The pressure distribution is then found to be 
C; ve b (* = 7) 4 
> - are Cos 
Cp:-dim. nV b? — | b —t 
c h-a-vVa-aa-2f| 
/ -In : (0) 
eV 1 —<c* t—¢c 


for -—1 S¢<1,b>1, -1<c <0. 
V?—1 


fort > 1,6 >1, —1<c <0. Inthe limiting case b = 1, 


Eq. (6) reduces to 


Cy ite, c 
C; = 4; —f tv) —¢ 


(6’) 


for —1<1< 1,0=1, -—1<c¢ <0 


By integration it was found that the lift coefficient 
based on the deflected area has its two-dimensional 
value. <A similar result was obtained for y = — 8 in 
reference 5. It follows easily from a_ reversibility 
theorem due to W. D. Hayes for the case when ly| > yw. 
For the case |/y| < yu, the explanation is to be found in 
another general theorem to be discussed elsewhere. 


TRIANGULAR PLAN FoRM: FREE EDGE SUBSONIC 


In the wing aileron combination shown in Fig. 4, 


upper and lower surfaces are no longer independent 














Fic. 4 


since the free edge FG is swept behind the Mach cone. 
This necessitates different methods of solution. The 
lifting case may no longer be reduced to the symmetrical 


case. However, the pressure distribution was easily 


obtained by the methods of conical flow. The results 


are :* 


wetes b—1)-V(b—o(1 +9) ‘ 


In 
“(1 + b) t—<¢ 
b lp — cy! ad 
—- 7) 
i+ 5V, —3u-u 
for—1<¢(<b,0<)<1,-1<c<0O. 


V1 — 2) 


. 9 
¢ P2-dim. ~ (= 


for —-1<ti<1b=1= —-e. 


In spite of the complicated formula for the pressure 
distribution, the center of pressure and the lift coeffi- 
cient have the simple forms 


(3b + 2c — 1)/4 (S) 


lop, = 


for0 <b<1,-18Sc¢ <0. 


* The method due to Evvard is also well adapted to this prob 
lem, especially when the edges of the control surfaces are curved 
Calculations carried out independently by H. S. Ribner, by this 


method, check the results given in the present report. 
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/ 
Ci/in = 4V (b — c)/(1 — 6) (9) 


forO0 SbS1,-1¢8c¢ <0. 


rhus, for a subsonic leading edge, the lift coefficient, 
based on the deflected area, is smaller than the two 
dimensional value. How it varies with the angles 3 
ind y is best seen from the accompanying graphs: 
Figs. 5-8. 
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Efficient Applications of Stringer Panel and 
Multicell Wing Construction 


GEORGE GERARD* 
New York Unwersity 


SUMMARY 


Under the structural loading conditions of transonic Mach 
Number flight, stringer panel type of wing construction and 
multicell construction were examined to determine the conditions 
under which each indicated efficient structural applications. 

From an analysis of structures of optimum design, stringer 
panel wings appear to have a wide range of efficient application. 
Only in the range of extremely high loadings and small thickness 
ratios does multicell construction indicate efficient applications. 

It is to be noted that the results obtained are strongly influ- 
enced by practical considerations involved in the manufacture of 


such wings. 


INTRODUCTION 


> EXTERNAL GEOMETRY of transonic aircraft wings 
is almost wholly dictated by aerodynamic require- 
ments of high Mach Number flight. It is necessary 
that compromises in aerodynamic design required to 
attain structural efficiency be kept to a minimum, 
and, therefore, the structure must be designed within 
limitations defined almost exclusively by aerodynamic 
considerations. 

Thus, design conditions are encountered wherein the 
thickness ratio of the wing is small and the bending 
moments large. Consequently, the cover must carry 
tremendous loading per inch. Also, it is essential that 
the cover be nonbuckling, and, while this requirement 
probably increases structural efficiency, it precludes 
construction that may cause serious shear lag effects. 

Under such conditions, it appears that in the pre- 
liminary design stages the designer has a choice of the 
type of construction (from a purely structural stand- 
point) which he may utilize for the wing. As shown in 
Fig. 1, it is possible to use skin-stringer panels for the 
covers separated by two webs which are more or less 
conventional design. It is also possible to use a multi- 
cell box in which a thick skin is stabilized by a series of 
webs. Presented with such a choice, it is of interest 
to investigate which type of construction has efficient 
structural applications from a weight standpoint. 


SYMBOLS 


Refer to Fig. 1: 

thickness of wing, in 
= empirical term, a measure of stringer panel efficiency 
E = modulus of elasticity, Ibs. per sq.in. 


u 


Received July 28, 1948 
* Instructor, Guggenheim School of Aeronautics, College of 


Engineering. 


K = instability coefficient, K = w*k/12(1 — »?) 

L’ = effective rib spacing (includes end fixity coefficient), in 
V; = loading per chordwise inch, lbs. per in. 

n = number of multicell bays (x + 1 = number of webs 
nby = structural chord, in 

t = thickness, in. 


= stress, lbs. per sq.in 
= = solidity 
y,¢, ¢ = variables defined in Eq. (17) 


cr. = critical 
w = web 
s = upper wing surface or compression cover 
t = lower wing surface or tension cover 
= — (over a symbol) denotes average or effective value 


ANALYSIS 


Assumptions 


(A) It is assumed that the modulus of elasticity, E, 
is constant, and no attempt is made to include correc- 
tions for inelastic buckling in the general case. 

(B) The wing is assumed to act under pure bending 
only, and, consequently, the structural weight required 
to carry shear loads is not considered. 

(C) The bending load is assumed to be resisted by the 
idealized boxes shown in Fig. 1, consisting of a “‘struc- 
tural chord”’ nb, and thickness by. 

(D) In order to have a reasonable criterion by which 
the upper and lower wing surfaces could be designed, 
it was decided that the lower surface would be designed 
to carry one-half the loading the upper surface carries 
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Fic. 1. Idealized box sections used for analysis 
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Fic. 2. Optimum curve for Y-stringer panels.” 

This would correspond approximately to an inverted 
flight condition, and, consequently, both surfaces 
would be designed from consideration of axial com- 
pressive loading. 


Optimum Panel Design 


If the wing utilizes skin-stringer panel construction 
as shown in Fig. 1, it is necessary to know, for the pur- 
poses of general analysis, the geometry of panels of 
minimum weight. 

The design of optimum panels was discussed by 
Zahorski,' who suggested that the optimum panel stress 
for a given compressive loading was obtained when in 
stability of the skin occurred simultaneously with pri 
mary and secondary instability of the stringers. For 
this condition, the optimum stress is given 


. , 1 
= CEN) (1) 


where c depends upon the crippling stress of the stringer 
and the panel geometry. 

The N.A.C.A.” has conducted an extensive program 
to determine optimum stringer panel proportions that 
would result in minimum weight for a given loading. 
The results of this program have shown that, of all 
shapes tested, use of a Y stringer panel results in mini- 
mum weight. These results are shown on Fig. 2. 

Also drawn on Fig. 2 is an extrapolation of the 
straight-line portion of the curve such as would prob- 
ably exist for an aluminum alloy with an extremely high 
proportional limit. It is interesting to note that the 
slope of this line is exactly one-half as predicted by Eq. 
(1),* although this line represents the envelope to a set 
of semiempirical curves obtained by varying the panel 
geometry. 

By taking the modulus of elasticity of the panels, 
which were constructed of extruded 75S-T Y-stringers 
and Alclad 75S-T sheet as EF = 10 X 10° Ibs. per sq.in., 
the equation of the optimum curve is 

@ = 1.152(EN,/L’)” (2) 

From Eqs. (1) and (2), c, which is the direct measure 


of the structural efficiency of the panel, has a value 


* The value of the exponent is obtained from the slope of the 


logarithmic plot. 
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for the Y-stringer panel of 
c = 1.152 3) 


By comparing values of this coefficient for various 
stringer shapes, a direct insight into the panel efficiency 
is obtained. 

From Zahorski’s analysis,' the coefficient ¢ has its 
highest value when the crippling stress and the radius of 
gyration of the stringer are highest and when the 
stringer spacing is close. The N.A.C.A.? investigation 
has substantiated this and has also shown that riveting 
of the stringer to the skin must be strong—4.e., large 
diameter rivets closely spaced. 

It is to be noted’ that, since the optimum curve is an 
envelope, its use ignores existence of only standard 
sheet gages and presupposes a large variety of extruded 
stiffener sizes. All of these considerations that lead to 
high structural efficiency present increasing difficulties 


in construction. 


Solidity of Stringer Panel Wing 

When considering the weight of a stringer panel wing 
section shown in Fig. 1, it is convenient to utilize the 
concept of “‘solidity.’’ Solidity is defined as the ratio 
of the area of the elements comprising the cross section 
to the total area of the cross section enclosed by the 
outer contours. It is evident that the solidity is di- 
rectly proportional to the weight of the structural 
material. 

In addition to the general assumptions of this analy- 
sis, the following assumptions pertain: (E) the weight 
of the ribs that are required to stabilize the cover is 
assumed negligible; (F) the centroidal depth of the 
stringer panel box, is assumed to be equal to the total 
depth, 6,.. 

The effective thickness of the upper stringer panel is 


Substituting for é from Eq. (1), 
, = (1/c)(NzL’/E)'” (5) 
According to assumption (D), 
Nz, = ('/2)N-z 
therefore 
I, = (1/c)(N,L'/2E)'” (6) 


The solidity of the box comprising the two covers and 
two webs shown in Fig. (1) is given by 


> = [(0, + nb, + Wyty|/nbdy (7 


Substituting Eqs. (5) and (6) for /., , and noting that, 
for a simply supported web in bending for which ¢-,. 


0, 
/ ” - ” . 
tie/ Dy = V12(1 — v*)&/25r°E 


the following equation for solidity is obtained 
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(1.707/c) (N,L’/E)' (2) 
b, 


/ 9 —_— pele ij 
— ”¢(N,/L'E)' ( =) (S) 
251? nb, 


Solidity of Multicell Wing 


nb, 


The method of analysis for determining the optimum 
number of webs to be used for a multicell wing has been 
presented in a previous paper.* The solidity of a multi- 


cell is given by 
> = [nbts + (n + 1)dutyo + nbgt,|/nbsdy (9) 


The thicknesses of the compression cover and webs 
are determined by the condition that the members are 


simply supported and buckle simultaneously. Thus, 
ty/ts = (7/5) (bw/ds) (10) 
and since 
Ns = Gat, = KjF2,°/0;* 
Thus 
te =VN,b,2/K,E (11) 


The thickness of the tension skin is determined from 


assumption (D). Thus, from Eq. (11) 


t,/ts = V'/2. = 0.794 (12) 

Combining Eqs. (10) to (12) with Eq. (9) as shown in 

reference 3, the expression for solidity of a multicell box 

that includes a lower surface designed to carry one-half 
the load of the upper surface becomes 


~ . BL Ne apf 1.794/nd,\ | (n+ 1)/ dy 
== | —— n + - (13) 
K,Enb, n b 5 nb, 


The optimum number of webs which will yield mini- 
mum solidity and consequently minimum weight is 


obtained from the condition 02/On = 0. Performing 
the differentiation leads to the result: 
n(4n + 1) = 5(1.794)(nb,/b,)? (14) 


In reference 3, the multicell was considered to have a 
lower surface designed to carry a tension load equal to 
the compressive load of the upper surface. Hence, the 
thickness of the tension skin was independent of the 
number of webs used, and for this case the optimum 
number of webs is given by 


n(4n + 1) = 5(nb,/b,)? (15) 


From Fig. 3, it is evident that use of Eq. (15) results 
in a smaller number of webs for a given structural thick- 
ness ratio, b,/nb,, than Eq. (14). Furthermore, use of 
the optimum number of webs given by Eq. (15) results 
in about a 4 per cent increase in solidity over that 
obtained by utilizing Eq. (14). Since the weight pen- 
alty is small for the large reduction in number of webs, 
it was thought desirable to use Eq. (15) in conjunction 
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Fic. 3. Optimum number of bays for a multicell box according 
to Eqs. (14) and (15) 
with Eq. (13) in this analysis, and this procedure was 


followed. 


EFFICIENT STRUCTURAL APPLICATIONS 


An examination of the solidity equations for the 
stringer panel and multicell box, Eqs. (8) and (13), 
indicates that the following variables are involved: 


Nz = cover loading per inch 
b,,/nb, = structural thickness ratio 
nb, = structural chord 

L’ = effective rib spacing 


(n, the number of bays, is specified by the structural 
thickness ratio for optimum conditions.) 


Consequently, for a given structural chord and effec- 
tive rib spacing, it is possible to construct a chart of 
solidity versus structural thickness ratio for various 
values of compressive loading. Such a chart was 
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Fic. 4. Solidity-thickness ratio chart for a specific case 


using optimum conditions. Intersection of solidity curves of 
stringer panel box and multicell determine boundary curve. 
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Fics. 5 and 6 


drawn in Fig. 4 for the specific case of mb, = 25 in., 
L’ = 24 in. 

In Fig. 4, the regions of efficient application of both 
types of construction are clearly evident. To the right 
of the intersection or boundary curve, a stringer panel 
box is more efficient, while, to the left of this curve, the 
multicell box is more efficient. 


Boundary Curve 


By equating Eqs. (8) and (13), the boundary curve 
can be directly determined without constructing a chart 
similar to Fig. 4. In doing so, the portion of the solid- 
ity contributed by the two webs in the stringer panel 
box may be neglected, since it was found by calculations 
that in the region of the intersection the webs contribute 
less than 5 per cent to solidity. Thus, equating Eqs. 
(S) and (13), 


(1.707 /c) (N,L’/E) (®) 
nb, by 


'y N, in| 1208 (abe) (2n + ” (2) | 
ess nN + 16 
K,Enb, nN b,, a nb 


Rearranging terms, 


(L")3 6 
NV, iD x 
(nh,)4 L.7O0TWRK, 
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Boundary curves that determine efficient structural applications of each type of construction 


Now, introducing the following 


y = (nb,)*/(L’)’ 
¢ = c/1.707 WK, 


‘2k Lae... BE Dro. ¥° 
p= n'[ + a+ (ee) 
n 5 nb, 


Eq. (17) can be written in the form 
N,/y7 = c8@*E (1s 

The ¢.term can be calculated from data of Fig. 3, 
since the optimum number of webs depends solely upon 
the structural thickness ratio. Thus, by letting c take 
on the value given in Eq. (3) and letting K, = 3.64, as 
for an efficiently designed simply supported plate of 
negligible taper,‘ and also letting H = 107 lbs. per sq.in. 
for aluminum alloys. Fig. 5 can be drawn. 

It is evident that the boundary curve is strongly in 
fluenced by the nature of c, since it appears as a sixth 
order term in Eq. (18). Thus, any departure from 
optimum panel dimensions such as use of standard 
sheet gages, wider spacing of stiffeners, and larger rivet 
pitch would result in a considerable shift in the bound 
ary curve. 

Similarly, use of less than the optimum number of 
webs would result in a shift in the boundary curve, since 
the ¢ term also appears as a sixth-power term. In 
reference 3, however, it was found that a reduction 


in the number of webs of two or three from the optimum 
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number resulted in an increase of only 5 per cent at low 
values of structural thickness ratios. 

It is interesting to note that the modulus of elasticity, 
E, enters into the equation for the boundary curve. 
Use of a steel alloy rather than an aluminum alloy 
would move the boundary curve to the left on Fig. 5, 
whereas use of a magnesium alloy would shift the curve 
tothe nght. It is not evident what effect exceeding the 
proportional limit would have on the boundary curve, 
particularly since the effective modulus is a function of 
stress rather than loading. 

In Eq. (18), the y term is a function of the structural 
chord and the effective rib spacing. To permit a better 
understanding of where the transition from stringer 
panel to multicell construction occurs, Fig. 6 was drawn 
for specific values of L’ and nb,. This figure indicates 
that, for optimum construction, the boundary curve 
vecurs at extremely low structural thickness ratios. 
Consequently, stringer panel construction can be uti- 
lized more efficiently than multicell construction in al- 
most the entire range of thickness ratios. 


DISCUSSION 


Che analysis conducted on structures of optimum 
design indicated that the transition from stringer panel 
to multicell construction occurs in possibly only a few 
highly specialized instances where the loading may be 
high or the thickness ratio may be low. The feasibility 
of constructing wings of the type specified by optimum 
considerations does not enter into the analysis, and, 
hence, constructional details may lead to conclusions 
that invalidate the result of the analysis in certain 
thickness ranges. Consequently, it is pertinent to con- 
sider what effect deviations from optimum construc- 
While the 
analysis involves certain ratios, it is necessary to con- 
sider the absolute dimensions of certain parameters in 
establishing limits. 


tion would have on the boundary curve. 


Stringer Depth 


An important consideration that would provide a 
cutoff to the boundary curve is that the stringer depth 
should probably not exceed one-quarter of the total 
depth of the structure. If the stringer were made con- 
siderably deeper, there is a definite loss in efficiency 
because the centroidal depth of the covers decreases, 
and also the stress distribution is such that the stress 
decreases from a maximum at the surface to zero at the 
center of the box. In addition, if the depth were too 
great, the stringers would touch. 


Riveting 

Rivet pitch and diameter have an appreciable effect 
on the strength of a stringer panel, particularly at high 
values of loading. According to reference 2, panels of 
optimum design, using N.A.C.A. countersinking meth- 
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ods, require an average rivet pitch equal to 6.5 times the 
skin thickness and a rivet diameter of 1.2 times the 
combined skin and stringer leg thicknesses. 

The rivet diameter requirement can be met, but the 
pitch requirement is such that the rivet heads are prac- 
tically touching. An investigation was made using the 
test data of reference 5 to determine the effect of in- 
creasing the rivet pitch to a value more likely to be used 
by the aircraft industry. For a rivet pitch equal to 15 
times the skin thickness, which appears to be a reason- 
able criterion based on industry practice, a reduction in 
failing stress, ¢, of approximately 15 per cent from the 
optimum value occurred for values of N,/L’ greater 
than 0.5 kips per in. 


Stringer Spacing 


Optimum conditions require that the skin and string- 
ers become unstable simultaneously necessitating close 
stringer spacing. If, for example, the centerline spac- 
ing of the Y stringers is increased to 4 in. from the opti- 
mum spacing, the data of reference 2 indicated a reduc 
tion of about 20 per cent in failing stress for a given 
N,/L’. 


Other Variables 


Other considerations that reduce the failing stress 
from the optimum value are use of standard sheet gages 
and only a limited variety of extruded Y stringer sizes. 
Since the Y stringer has the highest values of c, [Eq. 
(3)| and was used in the analysis, use of a Z stringer 
(which has the advantages that it can be formed, re 
quires fewer rows of rivets, and also results in better 
aerodynamic fairing) would result in a lower value of ¢ 
and a consequent reduction in failing stress. 


Multicell Webs 


From a manufacturing standpoint, it appears difficult 
to use more than three bays in a multicell wing. Con 
sequently, this factor tends to compensate at low thick- 
ness ratios for the above factors which tend to de- 
crease the efficiency of stringer panels. The effect of 
using only three bays can be estimated from refer- 
ence 3. 

‘It is evident that the range of variables that influence 
conditions other than optimum is so great as to pre- 
clude a general investigation for other than optimum 
conditions. For such cases, it is necessary to confine 
the study to specific cases. This can be done by con- 
structing a chart similar to Fig. 4 for the types of 
stringer panels and multicell construction considered 
feasible. 


CONCLUSIONS 


Based on optimum design of stringer panels and 
multicell wings and the assumptions of the analysis, it 
appears that stringer panel wings have a wide range of 
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efficient structural application. 
extremely high loadings and small thickness ratios does 


Only in the region of 


multicell construction indicate efficient application. 
However, the practical considerations involved in the 
manufacture of such wings influence strongly the bound- 
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ary between the two regions of efficient application. 
Consequently, an analysis such as presented herein 
should be made for the types of construction considered 
feasible to decide the type of structure to be used in 


specific cases. 
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The Boundary Layer of Yawed Infinite 


Wings 


J. M. WILD* 
Cornell Unwersity 


SUMMARY 


Au approximate solution to the boundary-layer equations for 
an infinite yawed cylinder is obtained by an extension of the 


Karm&n-Pohlhausen method. From this solution the direction 


of the flow in the boundary layer and the displacement and 
momentum thicknesses of the spanwise flow component may be 
calculated. 

To illustrate this method, the boundary layer over a yawed 
infinite wing at a lift coefficient just below the stall is calculated, 
and typical results are presented graphically. A strong span- 
wise flow is found which is in qualitative agreement with experi- 


mental observations. 


INTRODUCTION 


5 hen UNSATISFACTORY AERODYNAMIC CHARACTERIS- 
TIcs of sweptback wings at high angles of attack, 
coupled with their increasing use on tailless and high- 
speed aircraft, have emphasized a significant discrep- 
ancy between experimentally obtained results and 
theoretical predictions of these characteristics when the 
theory is based upon the frictionless-fluid approxima- 
tion. Accumulating experimental evidence—e.g., tuft 
studies indicating large spanwise outflows in the bound- 
ary layers of these wings—points to the strong possi- 
bility that viscous effects can no longer be safely ne- 
glected in sweptback or yawed finite-wing theories. To 
explore this possibility, a program has been laid out to 
investigate the three-dimensional boundary layers of 
finite wings and to determine the effects of such flows 
upon the wing characteristics. Sears,’ as a first step 
in this program, has calculated the cases of the laminar 
boundary layer of the yawed infinite flat plate and of a 
certain yawed infinite cylinder, both at zero lift. 

In continuation of this program, the laminar bound- 
ary layer of an infinite wing just below the stall is con- 
sidered here. It will be shown that an approximate solu- 
tion to this problem may be obtained by means of a 
simple extension of the Karman-Pohlhausen method. 


THE BOUNDARY-LAYER EQUATIONS 


It has been shown by Prandtl,? and independently by 
various other authors,’ * 4 that the Navier-Stokes 
equations (in the boundary-layer approximation) for a 
uniform, yawed cylinder infinite in extent in a steady, 
incompressible flow reduce to the form 

Presented at the Aerodynamics Session, Sixteenth Annual 


Meeting, I.A.S., New York, January 26-29, 1948. 
* Associate Professor of Aeronautical Engineering. 


Ou Ou 1 Op O*u mR 
u +w-=- - + » — (1) 
Ox Oz p Ox Oz? 
Op/dz = 0 (2) 
v v 07u 
+ w = yp (3) 
° On Oz Os? 


where x, y, 2 are curvilinear coordinates defined as in 
Fig. 1, » is the static pressure, p is the density, v is the 
kinematic viscosity, and u, v, ware the velocity compo- 
nents in the directions x, y, 3, respectively. In particu- 
lar, these authors point out that Eqs. (1) and (2) are 
independent of the spanwise velocity component and 
that these equations, along with the equation of con- 
tinuity, 

(Ou/Ox) + (Ow/Ozs) = 0 (4) 


constitute the usual two-dimensional problem of plane 
flow about the same cylinder at zero yaw in a free stream 
of velocity U’, the velocity component in the x-direc- 
tion. 

Approximate solutions to the plane problem were 
provided by Pohlhausen,’ who, at the suggestion of 
von Karman, expressed the boundary-layer velocity 
profile in the form of a quartic 


soe aGg)ea Gy ra)'s 
<a oe wes 
z 4 
a. (=), 0<2 
Or 


where u;, is the potential velocity outside the boundary 
layer, 6, is the boundary-layer thickness, and the A,’s are 
functions of x alone, determined to satisfy certain 
boundary conditions as well as Eq. (2). The boundary- 
layer thickness 6, is then calculated to satisfy Eqs. (1) 


6, (5) 


lA 








Sketch showing coordinate system used. 


Fic. 1. 
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and (4). This solution is well known and has been ap 
plied to a wide variety of problems with satisfactory 
results. 

Pohlhausen’s method may be extended to the de- 
termination of the boundary-layer velocity component 
in the y-direction on an infinite yawed wing i.e., to 
satisfy Eq. (3). To retain the same degree of approxi- 
mation as Pohlhausen, it will be assumed that the pro- 
file of the velocity component v7 can be expressed as a 


+ By (=) 
Oy 


quartic of the form 


~~ so +B (2) 
ay 8 by 


where v, is the local potential velocity outside of the 
boundary layer (in the case of the infinite wing, a con- 
stant equal to the y-component I’ of the free-stream 
velocity), and the B,’s are functions of x. The ‘“‘bound 
ary-layer thickness’ of the spanwise velocity compo 
nent, dy, is the distance over which v increases from zero 
at the airfoil surface to its potential-flow value, just as 
6, represents the corresponding distance for the com 
ponent w. For brevity, 6, will be referred to as the 
thickness of the “w-layer’’ and 6, as that of the ‘‘v 
10 


layer. 
¢ 6 
The functions Bo, B,, are now required to satisfy 


the following boundary conditions: 


v = 0 when z = 0 (Prandtl no-slip condition ) 

v” = 0 whens = 0 (ef. reference 6) 

” = v, when z = 6,)(To assure a smooth transi- ) (7) 
v’ = O when z tion to the potential ve- 


= §, 
v” = 0 when z = 5, | locity 


where the primes denote differentiation with respect to 
s. Under these conditions Eq. (6) becomes 


--2()-2() + (9 
vy by by by 


In this equation it is of interest to note that the non- 


us 8s 6, S 


dimensional profile of the v-layer is constant and, in 
particular, that separation profiles cannot occur (i.e., 
v’ # Oforz =0). Thisis consistent with the fact that 
for the infinite wing there are no spanwise pressure 
gradients. 

By integration of Eq. (3) with respect to z, or from 
momentum considerations, the momentum form of the 
v-layer equation is 


* To be more exact, of course, there exists only one boundary 
layer and a corresponding boundary-layer thickness over any 
body—namely, the layer wherein the resultant velocity vector 
deviates substantially from potential-flow vector 
nomenclature is used for convenience. The true boundary-layer 
thickness at any point of the surface is equal to the greater of 


The present 


bz and 6, 
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v : 
v Q)T 
Oz g=f) 


d ff 
u(v — v)dz 
dx 0 


If the values of w and v from Eys. (5) and (8) are in 
troduced into Eq. (9) and the integration carried out, 
an ordinary differential equation for 6,/6, results, the 
form of which depends upon the relative magnitudes 
then for 5, > 64 


of 6, and 6,. Let & denote 6,/6,; 


(€ Ss 1), 


and for 6, < 6,(€ > 1), 


ha it, 6. uy 
2 " a | T ¥Y 
, e , 
‘ r Uy Or uy 
¢ - (10Ob 


By u,’) 


where the prime denotes differentiation with respect to 


ry, and 
d u,'6,"/v, following Pohlhausen } 
A, B, C\ are functions of 4; B,, and &, as 
a, 38, yf — given in the appendix 


Ihe terms (11) are to be evaluated from the numeri 
cal results of the two-dimensional von Karman-Pohl 
hausen solution so that, in general, an analytical solu 
tion for £ is not possible. It can be seen, from Eqs 
(10a) and (10b), moreover, that at the forward stagna- 
tion line (x = 0) wm = 0 and the denominator vanishes. 
Thus, é’ is infinite at the stagnation line unless the 
numerator also vanishes. Rejecting an infinite value 
of &’ as physically impossible, one obtains from the ap- 
propriate Eqs. (10a) and (10b) a quintic equation giving 
the possible values of at the stagnation line. A nu 
merical solution indicates that the only physically ac- 
ceptable value is £ = 0.766. Thus, the initial condi- 
tion for the solution of Eqs. (10) is — = 0.766 at x = 0. 
Since this is a singular point, the best method for obtain 
ing a solution appears to be that of isoclines. An analo- 
gous situation exists in the two-dimensional problem 
and is discussed in various references.® ° 

From this solution, the v-layer profile is at once ob- 
tained, and such auxiliary results as the direction of 
flow in the boundary layer and the displacement and 
momentum thicknesses of the v-layer are readily calcu- 
lated. 
cated by the fact that a choice between Eqs. (10a) and 
The upper and lower parts of 


The calculations are not essentially compli- 


(10b) must be made. 
the isocline diagram are simply calculated from the 
respective equations, so that the shift is always made 
when the solution crosses the line & = 1. 


NUMERICAL RESULTS 


The configuration chosen to illustrate the above 
method is an infinite elliptic cylinder, with the ratio of 


t Cf. reference 7 
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Fic. 2. Lift curve of elliptical cylinder investigated. 


the inajor to the minor axes of 6: 1, at an angle of attack 
of 7° relative to the U-velocity component. This cor- 
responds to the maximum lift coefficient (shown in 
Fig. 2) as calculated by Howarth.* In carrying out 
this lift calculation, Howarth obtained the Pohlhausen 
solution for the laminar boundary layer over the un- 
yawed cylinder. These data were used in the subse- 
quent calculations and materially reduced the labor 
involved in obtaining the three-dimensional solution. 
In general, it would be necessary to obtain the two- 
dimensional solution before proceeding to the three- 
dimensional case. 

It was convenient to use the eccentric angle 7 (de- 
fined in Fig. 2) as the independent variable in carrying 
out the integration of Eqs. (10a) and (10b), rather than 
x. The integration of this equation by the isocline 
method is shown in Fig. 3, the solution being indicated 
by the heavy solid line. The dashed lines are two 
“neighboring’’ solutions constructed to check the ac- 
curacy of the isocline procedure but not satisfying the 
initial condition of finite —’ at x = 0. 

Typical profiles of the velocity components u and v 
corresponding to this solution are shown in Fig. 4. 
In Fig. 5 the ‘limiting streamline” at z = 0 is drawn on 
a projected plan view of a wing with 45° sweep, in 
comparison with the corresponding potential-flow 
streamline. Since the stagnation point for plane flow 
about the cylinder at this incidence is on the lower sur- 
face, the flow at this location is not shown in this figure. 
The short-line segments on this figure have been drawn 
in the “‘limiting’’ direction of the flow at the surface 
1 ¢€ 

dy/dx = lim (v/u) 
:—>0 

rhus, they may be imagined as representing yarn 

tufts attached to the surface of the airfoil to indicate 


the flow direction. For other angles of sweepback, the 
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flow direction may be obtained by multiplying these 
slopes by the tangent of this angle. 

In practice, the tufts undoubtedly assume a direc- 
tion corresponding to some intermediate streamline or 
streamlines. From the available data it appears that 
these intermediate streamlines lie between the limiting 
cases shown, but the w-velocity component would have 
to be determined in order to draw them accurately. 

The ratio of the displacement thicknesses 


bx 
as | [1 — (u/u)]dz 
Or 0 ; (12)° 


is shown in Fig. 6. Similarly, the ratio of the momen- 


tum thicknesses 


bx 
J (u/uy)(1 — (u/m) |ds 
6,, 0 . ‘ 
= - (13) 


6, oy 
ss f v/v)[1 — (v/m)]dz 


is shown in Fig. 7. 


CONCLUDING REMARKS 


A laminar boundary layer such as was assumed in the 
numerical calculations above could, in practice, exist 
only at extremely low values of the Reynolds Number 
because of the inherent instability of such layers. The 
results obtained here, nevertheless, appear to be in at 
least qualitative agreement with tuft observations on 
actual sweptback wings, since a strong outward flow in 
the boundary layer is obtained. This outward flow is 
more pronounced than that on the symmetrical cylin- 
der investigated by Sears,' presumably because of the 
stronger chordwise pressure gradients that are present. 
It would appear from this, at least qualitatively, that 
the outward flow is increased with increasing lift. This 
is also in agreement with tuft observations. 

It should be noted, however, that this spanwise com 
ponent of flow does not have any effect upon such 
aerodynamic characteristics as the lift, transition, or 
separation for the infinite wing. In fact, the basic 
Eqs. (1)-(4) are such as to preclude the possibility of 
any such effects occurring.* It appears that these ef- 
fects can arise only on finite wings. 

The present extension of Pohlhausen’s procedure to 
the yawed wing appears suitable for estimation of the 
laminar boundary-layer characteristics in certain areas 
of sweptback wings. In view of the analogy between 
the v-layer and the temperature layer near a heated 
surface,'! the same procedure should be useful in analo- 
gous heat-transfer calculations for either straight or 
swept wings. 





* A more detailed discussion of this point is given in references 


l and 4 
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Typical spanwise < 


and chordwise velocity components in the boundary layer at varying chordwise locations along the 
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Plan view of the potential flow and limiting streamline 
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APPENDIX 


in Eqs. (10a) and (10b) may be 
evaluated from the following: 


The coefficients 


A = 0.3§-! + D, — Dot + D3t? — Dyt4 
B = —0.3¢-? — D. + 3D,;&? — 4D,¢° 
C = D,’ — D,/§ + D3’ — Dy’ 
where 
4 
D, = > [A,/G+ 1) -1 
t=1 
1 
D, = >> [2A,/(i + 2)] - 1 
t=] 
4 
D; = > [2A:/(i + 4)] — (1/2) 
$=] 
4 
Dy = > [Ai/(i + 5)] — (1/5) 
t=] 
and 
a = 0.06667A;, + 0.02381A.¢-! + 0.01072A;¢-? + 
0.00556A ,E-8 
8 = —(0.13334A,E— + 0.07143A.-? + 0.04288A3£-8 
+ 0.02780A,é-*) 
y = 0.06667A,’ + 0.02381A.’—-! + 0.01072A;’E-? + 


0.00556A 4’&~* 


The primes denote differentiation with respect to x. 
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The B,’s of Eq. (6) have been introduced into these ex- 
pressions numerically as indicated by Eq. (8). 
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Problems and Progress in Low-Pressures 
Research 


E. D. KANE? anpb R. G. FOLSOM? 
Unwersity of California at Berkeley 


ABSTRACT 


The nature of the theoretical and experimental problems 
associated with the field of rarefied gas dynamics is discussed. 
The types of laboratory equipment judged to be most suitable 
for studying these problems are low-pressure supersonic wind 
tunnels, molecular beams, and rotating cylinders. Experimental 
results with a model wind tunnel and typical supersonic nozzles 
illustrate the existence of the expected large viscous effects and 
possible slip at the boundary. The techniques used and under 
development for pressure measurement and flow visualization 
are described. 


INTRODUCTION 


bien MECHANICS OF FLOW of a gas at low pressures, 
or rarefied gas dynamics, has reached the stage of 
engineering importance. The term “low pressures’ 
is understood to mean low absolute pressures of | in. 
Hg or less. Subsonic gas flow in pipes under low- 
pressure conditions prevails in vacuum systems used 
in the preparation of vitamins, drugs, or in similar 
distillation processes. There is a fair body of theory 
and experimental information available for this type 
There are virtually no experimental data or 


adequate theoretical treatments available for super- 


of flow. 


sonic flow in a rarefied gas (sometimes called ‘‘super 
aerodynamics’), and it is the investigation of this 
field which is the primary purpose of the research 
program to be discussed. For the aeronautical engi- 
neer or the aerodynamicist, a knowledge of fluid me 
chanics as applied to the supersonic flow of a gas at 
low pressures is necessary to understand the problem 
of flight at extreme altitudes. 

Applications of the available theory of gas dynamics 
to low-pressure phenomena are limited because the 
fluid cannot be described in terms of a continuum but 
must be discussed on the basis of the motion of indi 
vidual malecules. A useful criterion in estimating the 
magnitude of low-pressure effects is the molecular 
mean free path, which can be defined approximately 
as a distance a molecule will travel, on the average, 
before it undergoes a collision with another molecule. 
When the pressure is high (or more properly speaking, 
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the density), the individual molecules of gas are close 
together, and the mean free path is correspondingly 
small. When the pressure is decreased, the number 
of molecules per unit volume decreases, and the average 
distance that an individual molecule will travel before 
impact increases. In order to fix the magnitude of the 
mean free path quantitatively, Table 1 shows some 
values of molecular mean free path at different altitudes 
and pressures. For example, at an altitude of 75 miles 
the molecular mean free path is about | ft. At this 
altitude, the atmosphere can no longer be considered 
a continuum for an object with a characteristic dimen- 
sion of 1 ft. If the characteristic dimension of the 
object was 1,000 ft., then the atmosphere might behave 
essentially like a continuum for that large object. 


Elementary physical reasoning suggests that the 


ratio of mean free path divided by characteristic 
dimension (1/D) provides one important parameter 
to describe the magnitude of low-pressure effects. .\s 
discussed by Tsien,’ this L/D ratio is proportional to 
the ratio Ma/Re, the constant of proportionality being 
roughly 1.5. The question of how large the L/D or 
Ma/Re number ratio must become before the methods 
of continuum fluid mechanics begin to require correction 
is not clear. Based on theoretical and experimental 
considerations for subsonic flow in a pipe, the ordinary 
gas dynamics theory requires correction when the 
L/D ratio becomes larger than about 0.001. Scham- 
berg’ estimates that if L/D is 0.005, the skin friction 
coefficient for a flat plate moving at a Mach Number 
of 3, decreases by 10 per cent from the value predicted 
by ordinary gas dynamics theory. 


As pointed out by Tsien,’ each phase of the non- 
continuum flow problem can be subdivided logically 
into three parts, depending on the magnitude of the 
pressure. For appreciable but small mean free paths, 
called the slip region, ordinary gas dynamics begins to 
fail. For mean free paths that are large compared 
to the body dimensions (L/D large), collisions of a 
molecule with a surface are much more frequent than 
collisions with other molecules. This region might 
be termed ‘‘fully developed’’ molecular flow. There is 
a transition region where the frequency of collisions 
with a wall and collisions with other molecules is of the 
same order of magnitude. Analytical treatment of 
these three regions may well require three kinds of 
theoretical development. Large gaps exist at present 





Cha 


in 
met 


a 
con 
flov 
ide 
bou 
are 
mol 
The 
nor 
At 
Rey 
dec 
for 
pre: 
bou 


whe 


The 
the 
niti 
be» 


Div 


For 


Eq. 
exp 
equ 
thi: 


close 
ingly 
nber 
Tage 
fore 
r the 
ome 
udes 
niles 
this 
ered 
nen- 

the 
lave 


the 
istic 
eter 

As 
1 to 
ing 
) or 
ods 
‘ion 
ital 
ary 
the 
Ln 
ion 
ber 
ted 


on- 
lly 
the 
hs, 

to 
‘ed 
Zz 
an 
ht 


ns 
he 
of 
of 





LOW-PRESSURES 


RESEARCH 47 


TABLE 1 


Molecular Mean Free Path for ‘Adopted Atmosphere” and for Room Temperatures 





——‘‘Adopted Atmosphere’’*——— 








—_—_—Pressure—— —— Altitude Temperature, Mean Free Mean Free Path 

Mm. Hg Lbs. per sq.ft. Miles Ft. X 10-5 “7. Path, Ft. at 70°F., Ft 

760 2,110 0 0 70 3.1 xX 3.1 X 107? 

180 500 6.6 0.35 —70 1.0 X 10-* 1.3 X 1078 

7.9 22 19 1.0 —70 2.5 X 10-5 3.0 X 1075 

0.57 1.6 32 1.7 170 5.0 X 1074 4.1 X 107' 

2.4 KX 107° 0.075 49 2.6 —30 7.5 X 107% 8.7 X 107° 
3.6 X 1074 0.99 * 10-3 76 4.0 220 1 0.66 
6.5 X 1076 1.8 xX 10 114 6.0 760 90 36 
5.0 X 10-° 1.4 X 10-¢ 152 8.0 1,300 1,500 460 


* Values from AA F Technical Report 5632, ‘‘A Design Manual for Determining the Thermal Characteristics of High Speed Aircraft"’ 


Chapt. III, 1947. 


in portions of the theory and particularly in experi- 
mental information. 
LAYER 


THE BOUNDARY 


. 

Shortly after 1900, Prandtl introduced the modern 
concept of flow of real fluids, pointing out that the 
flow of a real fluid of small viscosity differs from an 
ideal fluid only in the region near a solid surface (the 
boundary layer). In this region, the frictional effects 
are appreciable, and the fluid undergoes a loss of 
momentum related to the frictional force at the surface. 
The general character of the boundary layer is indicated 
normally by the magnitude of Reynolds Number. 
At low pressures the boundary layer is laminar, since 
Reynolds becomes because of the 
decrease in density. Solutions have been developed 
for laminar flow over a smooth flat plate with zero 
pressure gradient in the of flow. The 
boundary-layer thickness 6 becomes 


Number small 


direction 


1 
6 = (constant) (xv/u) ” (1) 


where 


’ = distance from leading edge of plate 
vy = kinematic viscosity of fluid 


u = free-stream velocity of fluid relative to plate 


The constant is estimated from 1.7 to 8, depending on 
the assumptions as to compressibility, the exact defi- 
nition of 6, etc. In order of magnitude, Eq. (1) can 


be written 


6 = SV xv/u (2) 


Dividing by a characteristic linear dimension D gives 


6/D = 5Rep~'* Vx/D (3) 

For 6/D = '/2, let x = xo in Eq. (3) and solve for x. 
%)/D = Rep X 107? (4) 
If D is the diameter of a nozzle and Rep = 100, 


Eq. (4) indicates that the boundary layer may be 
expected to fill the nozzle when the nozzle length 
equals the nozzle diameter. Reynolds Numbers of 
this order are possible in small nozzles at low pressures, 


and it is apparent that viscous effects become im 
portant to flow geometry at low pressures. 

When mean free paths are large, experiments show 
that another effect occurs which modifies Eqs. (1) to 
(4). One of the boundary conditions used in deriving 
Eq. (1) was that the fluid velocity at a solid surface is 
zero. Low-pressure flow experiments in pipes indicate 
that a finite “‘slip’’ velocity occurs at a solid boundary, 
resulting in a change of the velocity distribution from 
the surface out to the free stream. The “‘slip”’ 
surface phenomenon, but, unlike the boundary layer, 
its character is determined by the nature of molecular 
interactions at the surface rather than by the viscous 
characteristics of the fluid. Information regarding 
reflection (specular or diffuse) and accommodation 
coefficients becomes of primary importance. A solu- 
tion for flow over a flat plate may be obtained under 
the assumption that the velocity distribution function 


is a 


remains the same as in laminar flow, except that a dis 
continuity occurs at the boundary. This solution 
gives an asymptotic expression for “boundary-layer”’ 
thickness as 


6/x ~ L/x for L/x > 1/Ma (5) 


where L = molécular mean free path, Ma = Mach 
Number and the other quantities were defined pre- 
Eq. (5) indicates that at low pressures the 
independent of 


viously. 
boundary-layer thickness 
Reynolds Number, a result quite different from Eq. (1). 

Present estimates of the magnitude of L/D for the 
“slip’’ region depend on experiments in pipes where 
no curvature, small pressure gradients, and no shocks 
Preliminary experiments with nozzles indicate 
Similar 


becomes 


exist. 
that ‘‘slip’’ may occur at smaller L/D values. 
results might be expected with airfoils when curvature 
and pressure gradients are appreciable. 

In supersonic flow the nature of the shock wave and 
shock patterns are important. The shock wave is 
thin at normal pressures, although elementary con 
siderations would indicate that the minimum thickness 
is probably several molecular mean free paths. As 
the density decreases, the shock-wave thickness must 
increase and the usual concept of the shock as a dis- 


continuity no longer holds. The interactions of a 
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shock region at low pressures with a laminar boundary 
layer that is partly subsonic and partly supersonic are 
complex. 


POSSIBLE TYPES OF EXPERIMENTAL EQUIPMENT 


With the problem of fluid flow at low pressures 
roughly defined, the question of suitable types of 
laboratory equipment for experimental investigations 
arises. In general, studies may be made on objects 
moved through stationary air or on stationary objects 
placed in an air stream. Special equipment will be 
necessary to study the details of the boundary phenom- 
ena. 


(1) Flight Testing 


Certain types of flow data may be obtained from 
recording or telemetering instruments carried by a 
rocket into the upper atmosphere as in the experiments 
being conducted with the V-2 at White Sands.* The 
highest altitudes reached have been in the order of 
115 miles. 


(2) Towing the Model 


Since large volumes at low pressure are achieved 
easily with modern vacuum pumps and design tech- 
niques, it is reasonable to consider sending a model 
through such an atmosphere. A firing chamber (of 


the type available at Aberdeen Proving Ground) 


could be adapted for this purpose. A variation of this 
idea is to attach a model to the end of a rotating arm 
that rotates at high speeds in a minimum-sized tank. 
The major difficulty with these schemes is instrumenta- 
tion. Drag forces at low pressures are small, and a 
fast recording of pressure readings is difficult at low 
pressures. 


(3) Low-Pressure Wind Tunnel 


A conventional supersonic wind tunnel with modi- 
fications designed to give low-pressure high-speed flow 
is an attractive possibility for studies in the ‘“‘slip”’ 
region. 
can be studied conveniently. The instrumentation 


The flow about different stationary models 


problems are simplified, and the operating conditions 
are under the control of the investigators. 


(4) Molecular Beam 


A molecular beam (of the type used for atomic and 
nuclear studies) provides a uniform stream of molecules 
suitable with 
Information on 


for studies of molecular interaction 
surfaces at extremely low pressures. 
reflection and accommodation coefficients 


obtained for controlled conditions at a given surface. 


may be 


(5) Rotating Cylinder 


Tsien' suggested that an extension of Millikan’s 
rotating cylinder experiments to high speeds would 
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provide fundamental information on “‘slip’ flows. 
Preliminary design studies have indicated that with 
suitable geometry the whole range of molecular flow 
phenomena, through “fully developed 


molecular flow’’ 


from 
could be obtained with this type of 


“slip” 
equipment. 


INSTRUMENTATION PROBLEMS 


In considering the use of experimental equipment 
to investigate low-pressure problems, the question of 
Static or 
impact pressures can be measured with an appropriate 


suitable instrumentation is important. 


type of low-pressure gage.’ In dealing with gas-flow 
problems, it is necessary that at least one other property 
of the gas be measured. Temperature is a possibility, 
but the difficulties in measuring local temperatures are 
great. 
more attractive possibility, particularly if the means 


Direct measurement of density seems to be a 


used to obtain quantitative density measurements 
can also be used to show visually the nature of the 
flow pattern. 


(1) Pressure Measurements 


If pressures at a wall tap or from an impact tube 
are to be obtained under low-pressure conditions, it 
is necessary to choose the dimensions of the system 
with care if impractically long response times are to 


be avoided.® Fig. 1 illustrates a typical gaging system 
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idealized as an orifice (wall tap), connecting line, and 
gage tube volume. The curves indicate the effect of 
geometry on the time required to reach equilibrium 
after the system pressure is changed. Response times 
in the order of hours are easily possible with improper 
designs. 

The meaning of an impact tube reading when the 
mean free path is of the same order as the tube diam- 
eter is a problem under investigation. A theoretical 
analysis has been developed and an experimental check 
will be made. 

Three basic types of pressure gages are being used 
on wind tunnel measurements. The McLeod gage‘ 
is the standard, and its use is limited to secondary 
gage calibration. It can be constructed to cover 
almost any low-pressure range, but considerable time 
is required to obtain a reading. A precision U-tube 
manometer, filled with butyl phthalate mixed with low 
vapor pressure oil to give a convenient specific gravity 
at 70°F., is used for pressures over 500 microns (0.5 
mm.) and up to about 6 mm. of Hg. The Pirani 
gage’ is useful over a pressure range from 1 to 500 
microns of Hg. In the Pirani gage, the temperature 
of a hot wire is a function of the pressure. Normally, 
two Pirani tubes are provided, one of which is sealed 
off at a definite pressure and placed in one leg of the 
bridge circuit, while the other is connected to the 
system. If the two tubes are maintained at the same 
temperature, the effect of atmospheric temperature 
variations is essentially canceled out. An indicating 
and recording system capable of giving the readings 
of 16 gages simultaneously has been developed. 

Probably the most difficult problem pf pressure 
measurement is the detection system of the molecular 
beam. With a background pressure of 5 X 10-* mm. 
Hg, it will be necessary to detect a pressure differential 
of 10-"'mm. Hg. The first development in this direc- 
tion uses ion gages (triode tubes in which pressure 
changes produce changes in plate current), one tube 
being used to probe the beam, while the other registers 
background pressure. Pressure differences of 5 X 
10-° mm. Hg have been measured with present equip- 
ment. 


(2) Flow Visualization 


For ordinary wind-tunnel work the shadowgraph, 
Schlieren, and interferometer techniques are used to 
visualize the flow. These light refraction methods are 
subject to decreased sensitivity as the fluid density 
decreases. At the lower pressures (less than 1 mm. 
of Hg abs.) these methods can no longer be employed 
in small wind tunnels. 

In attempting to develop new methods of flow visuali- 
zation, consideration has been given to the property 
of molecules for absorption of radiation of given wave 
lengths or emission of radiation when properly excited.® 
Typical possibilities are exhibited by oxygen, which 
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has a strong absorption band at a wave length of about 
1,470A.; by mercury vapor with a strong absorption 
band at 2,537A.; and by ozone, which has its maximum 
absorption at about 2,537A. Preliminary investiga- 
tions have shown that flow visualization is possible 
for a limited range of low pressures with any of these 


gases. In general, three problems are involved: 


(1) The generation of monochromatic light of the 
proper wave length. 

(2) The production and introduction of the gas or 
vapor. 

(3) Detection of the radiation passing through and 
absorbed partially by the gas. 


Preliminary work with the mercury vapor absorption 
flow visualization process indicates expected difficulties 
connected with low vapor pressure of mercury and the 
proper interpretation of measurements because of the 
reactions of the relatively heavy mercury molecule. 
The quantities of ozone required for reasonable sensi- 
tivity with an ozone absorption method are of such 
magnitude that ozone generation and storage become 
a major problem. Ozone in liquid form is unstable 
and highly explosive. The oxygen absorption flow 
visualization method appears to offer the best possi- 
bilities. 

The oxygen absorption method has been studied in 
detail. A hot cathode discharge tube with a Xenon 
atmosphere has been developed for the production 
of light of the required 1,470A. wave length. Details 
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of an experimental lamp are shown in Fig. 2. The 
cathodes are fed from a 220-volt, 60-cycle supply. 
Starting of the tube is accomplished by means of a 
radio-frequency (27.4 me.) field at the capillary (ap- 
proximately !/, in. I.D.) of the tube. 

As indicated in Fig. 2, the capillary section is viewed 
from the end through a fluorite window. This par- 
ticular tube was provided with a chamber around the 
window in which the oxygen pressure could be varied 
independently of the tube pressure. Radiation passing 
through the fluorite window excited fluorescent ma- 
terial coated on a second window in the chamber. If 
the oxygen pressure in the chamber was raised, some 
of the radiation was absorbed, and the fluorescent 
intensity decreased. By measuring this intensity with 
a photomultiplier tube as a function of oxygen pressure 
in the chamber, the absorption cross section of oxygen 
for the lamp radiation was determined. This indirect 
procedure was followed because of unavailability of 
a vacuum spectrograph. The experiment showed that 
good intensities of 1,470A. light were emitted by the tube. 

In applying the method to the wind tunnel, the 
1,470A. source would be used with a suitable fluorite 
optical system to give parallel light at the test section. 
Light traversing the test section would be absorbed by 
an amount depending on the path length and density, 
giving a nonuniform light distribution at exit from the 
test section. The changes in intensity over the test 
section would then be detected by photographic films, 
observations of a fluorescent screen, or equivalent 
means. To indicate the magnitude of the effect, for 
a 2 to 1 density ratio, a change in 1,470A. intensity 
of 50 per cent is expected at 100 microns of mercury 
if the gas is 100 per cent oxygen. If air is used, the 
intensity change will be about 10 per cent. 

When properly excited, nitrogen produces a radiation 
that persists after the excitation has been removed. 
This afterglow phenomenon, being explored by Prof. 
J. Kaplan, may provide a practical method for flow 
visualization at low pressures. 


Low-PRESSURE SUPERSONIC WIND TUNNELS 


The development of wind tunnels at Berkeley for 
low-pressure supersonic flow considers only the con- 
tinuous operating, nonreturn, open, or closed throat 
type of tunnel. Three of tunnels are under 
investigation, the smaller units being limited in use- 
fulness but more economical to operate. Preliminary 
investigations with the smaller units have proved 
necessary in determining the performance and geometry 
of the largest tunnel. 

The smallest tunnel’ has a 1- by */s-in. rectangular 
test section operating at a pressure of 1 in. Hg abs. 
at M = 2.4. A single-stage steam ejector exhausting 
to the atmosphere provides the driving energy. This 
tunnel will not operate at pressures low enough to show 
molecular flow effects. Preliminary work has been 


sizes 
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completed, and the unit now provides equipment suit- 
able for many student tests. 

The largest tunnel now under construction will have 
a test section of about 150 sq.in. in area and will be 
suitable for tests on models with a characteristic di- 


mension of about | in. 


The pumping system is based 
on a five-stage steam ejector’ selected on the basis of 
economic operation and flexibility of performance. 

The medium tunnel was built with a I- by 1-in. test 
section to permit experimentation on nozzle and diffuser 
design on an economical scale and to provide a suitable 
means for developing instrumentation for low-pressure 
conditions. The general arrangement of the tunnel 
is shown in Fig. 3. Tracing the flow cycle of gas 
through the tunnel, atmospheric air is passed through 
a dryer, flow rate meters, a needle valve, and into a 
settling chamber or reservoir where stagnation pressure 
and temperature can be measured. The gas then 
passes through a supersonic nozzle to the test section; 
then through a diffuser into a large ballast tank, the 
purpose of which is to damp out high-frequency fluctua- 
tions in the pumping system performance; and then 
into the first-stage vacuum pump. This pump raises 
the gas pressure to the range of 200 to 3,000 microns 
of Hg, depending on operating conditions. The gas 
is raised to atmospheric pressure and discharged by a 
second pump of the mechanical positive displacement 
type. 
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No. 2 wind tunnel (schematic sketch combined with 
photograph). 
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It will be noted from Fig. 3 that a high vacuum 
pump (capable of producing pressures of about 10~° 
mm. of Hg abs.) is connected to the ballast tank. This 
pump is useful for calibrating the pressure gages with- 
out removing them from the tunnel. The calibration 
standard is a McLeod gage. 

To start the tunnel, the inlet needle valve and down- 
stream valve are closed, and the mechanical pump is 
connected directly to the wind tunnel. When no leaks 
exist this pump will bring the pressure in the tunnel 
from atmospheric down to about 0.75 mm. Hg in 3 
min. The pump is then closed off from the tunnel, 
the downstream valve is opened, and the fore pump in 
series with the mechanical pump brings the pressure 
down to less than 5 microns of Hg in about 2 min. 
Then the inlet valve is opened sufficiently to give the 
desired rate of gas flow through the tunnel. Control 
of tunnel operating conditions is obtainable by adjust- 
ing one or more of the following variables: (a) inlet 
valve position, (b) downstream valve position, and 
(c) pump performance. 

The volumetric and weight rate of flow obtainable 
with the first-stage pump at a given intake pressure 
can be varied by changing the power input to the 
pump. This high vacuum pump is a commercial oil 
jet unit, using a low vapor pressure oil that is vaporized 
in a boiler and expanded through a supersonic nozzle 
to give a pumping action analogous to that of a steam 
jet ejector or water jet pump. Approximately constant 
volumetric flow rate can be obtained with this unit 
over a pressure range from about 70 to 300 microns 


of Hg. 


PRELIMINARY EXPERIMENTAL RESULTS 


1) Diffuser Performance 

The viscous effects at low pressures make the 
problem of diffuser design for pressure recovery diffi- 
cult. It is probable that large pressure ratios, of the 
magnitude obtainable in supersonic tunnels operating 
at normal pressures, cannot be achieved in a wind 
tunnel for low-pressure testing. Tests were run with 
the diffuser shown in Fig. 4 in position in the tunnel. 
The position of the diffuser along the axis of the tunnel 
was adjustable in the range from about !/s to 2 in. 
from the end of the nozzle. With this particular diffu- 
ser the ratio of ballast tank pressure to nozzle static 
pressure was about 1.3 with the diffuser in place as 
compared to 1.1 with no diffuser. Because of the 
small gain in tunnel performance achieved in this 
diffuser, it has not been used for regular test work. 
Development of a possibly better diffuser has not been 
considered of major importance at this time. 


(2) | Wind-Tunnel Performance 


The performance of the wind tunnel with a par- 
ticular nozzle is shown in Fig. 5. These data were 
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Fic. 4. Photograph of nozzle in No. 2 wind tunnel. 
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obtained by setting the inlet valve to give a constant 
reservoir or supply pressure to the nozzle of 5 mm. of 
Hg. The downstream valve was used to throttle the 
flow and control the test chamber and ballast tank 
pressures. Impact tube pressures were read at the 
center of the exit plane of the nozzle and were com- 
bined with static pressures read at the nozzle wall to 
calculate the apparent Mach number, J. 

The first effect noticed when operating the tunnel 
was that, in closing the downstream valve, the pressure 
in the test section as measured at the nozzle wall tap 
immediately began to increase. This phenomenon 
was contrary to what would be expected in a super- 
sonic wind tunnel operating at normal pressures, 
where the effect of downstream changes is not apparent 
at the test section until a shock pattern has been moved 
forward into that section. The tentative explanation 
is that the thick boundary layers existing in the tunnel 
permit transmission of downstream effects to the up- 
stream portions of the nozzle through fluid which is at 
subsonic velocities. 

Using the one-dimensional nonviscous compressible 
flow theory, a constant Mach Number = 4 would be 
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expected with this nozzle. From Fig. 5, the maximum 
value of Mo for the given reservoir pressure (with the 
downstream valve open) was about 2.55. It is charac- 
teristic of the operation of this low-density wind tunnel 
that the Mach Number obtainable depends on the 
effective sectional area of the nozzle as influenced by 
the boundary-layer thickness. 


(3) Nozzle Performance 


Using an impact tube (0.027 in. I.D. and 0.047 in. 
0.D.), pressure distributions were obtained at a plane 
\/, in. upstream from the exit plane of a rectangular 
nozzle (designed in the diverging portion by the method 
of characteristics for JJ = 4 with no correction for 
boundary-layer thickness). The results are expressed 
in terms of an apparent Mach Number determined by 
the usual Rayleigh formula, using the static pressure 
measured at a pressure tap at the plane of the impact 
tube in a diverging wall of the nozzle. In all tests the 
diffuser was removed, and the control of the pressure 
was by means of the tunnel inlet valve. Theoretical 
and experimental work indicate the apparent Mach 
Number as used is fictitious because of static pressure 
variation in the test section and the influence of large 
molecular mean free path. The experimental data are 
presented in this form for convenience. Additional 
work is in progress to determine the magnitude of the 
corrections that should be applied. 

A typical impact pressure distribution along the 
vertical axis (between diverging walls) is shown in 
Fig. 6. The characteristic thick boundary layers 
and lack of slip at the walls are noticeable. All plotted 
data are uncorrected, and the apparent shift in values 
at the centerline axis is due to experimental runs in 
each half of the nozzle being made on different days 
with slightly different adjustments of the tunnel. 
Fig. 7 shows similar data for the horizontal axis at a 
series of static pressures. These data indicate the 
possibility of slip at the parallel walls of the nozzle. 
For small openings of the inlet valve, the pressure is 
low and the Mach Number is correspondingly low, and 
the distribution is characteristic of what might be 
expected for laminar flow at subsonic speeds. As the 
inlet valve opening is increased, the pressure rises, the 
Mach Number increases, and a uniform velocity 
section appears at the center of the nozzle, becoming 
about 25 per cent of the cross section at an indicated 
static pressure of 200 microns. The presence of a thick 
boundary layer and large viscous effects is demon- 
strated. 

The major problem in designing a supersonic wind 
tunnel for low-pressure studies is to determine the 
type of nozzle that will give optimum velocity dis- 
tribution. Shortening the nozzle length in the flow 
direction in an attempt to reduce boundary-layer 
build-up is one obvious approach that is now being 
examined both experimentally and theoretically. The 
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Fic. 8. Molecular beam (schematic sketch combined with 
photograph). 


boundary-layer thicknesses indicated (Figs. 6 and 7) 
may be due partially to separation effects. Before 
constructing the nozzle with the abrupt converging 
section shown in Fig. 4, unpublished tests were run 
with various geometries. Oil was used at velocities 
to give the Reynolds Number range encountered in 
the wind tunnel, and velocity distributions were ob- 
tained at each nozzle throat. These tests with oil 
showed that the abrupt approach section used in the 
wind tunnel supersonic nozzle gave a uniform velocity 
distribution at the throat with minimum boundary- 
layer thickness. Additional evidence of satisfactory 
operation was available from performance of similar 
sharp entrance nozzles used in the Kochel wind tunnel. 
It is possible that shocks occurring in the supersonic 
portion of the nozzle would cause thickening of the 
boundary layer observed near the nozzle exit. 


MOLECULAR BEAM 


Molecular beam equipment (Fig. 8) is being tested 
at Berkeley. The source of high-speed molecules is a 
small furnace with an opening in one end through 
which molecules can escape into an evacuated chamber. 
This ‘“‘beam’’ of molecules is collimated with slits to 
give a cone of molecules traveling with a range of speeds. 
A velocity selector can be used to pass only molecules 
having speeds within a predetermined range. The 
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instrument is under development with the ultimate 
objective of determining the nature of the interaction 
of molecules with surfaces. Molecular (N2) speeds 
of about 1,600 ft. per sec. are used currently in equip- 
ment designed for a maximum of about 3,500 ft. per 
sec. at high furnace temperatures. 
, ‘“‘model.”’ 


The beam is about 
is in. diameter at the 
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Airfoils in Slightly Supersonic Flow 


* 


HSUE-SHEN TSIEN? anp JUDSON R. BARON? 
Massachusetts Institute of Technology 


SUMMARY 


An investigation is made to determine the performance of 
simple thin airfoils in the slightly supersonic flow region with the 
aid of the nonlinear transonic theory first developed by von 
K4rm4n.! Expressions for the pressure coefficient across an 
oblique shock and a Prandtl-Meyer expansion are developed in 
terms of a transonic similarity parameter. Aerodynamic coef- 
ficients are calculated in similarity form for the flat plate and 
asymmetric wedge airfoils, and curves are plotted. Sample 
curves for a flat plate and a specific asymmetric wedge are plotted 
on the usual coordinate grid of C;, Cz, and Cres, versus angle of 
attack and C; versus Mach Number to illustrate the apparent 
features of nonlinear flow. 


INTRODUCTION 


gegen THEODORE VON KARMAN! has derived a 
set of similarity rules for flows around slender 
bodies in a stream having close to sonic velocity. He 
called this set of rules the transonic similarity rules. 
Essentially, the results are as follows: If 6 is the thick- 
ness ratio of the body, a the angle of attack of the air- 
foil, C, Ca, Cm the sectional lift, drag, and moment 
coefficients, respectively, then 


Ci(a/s) = (6/T'")L(K) (1) 

Ca(a/6) = (6°/T'")D(K) (2) 
and 

Cn(@/6) = (6/0) M(K) (3) 


In these equations, K is the similarity parameter de- 


fined as 


K = (1/2)[(M°? — 1)/(T8)”] (4) 
where M° is the free stream Mach Number and I isa 


parameter of the property of the fluid—i.e., 


r= (y+ 1)/2 (5) 
where vy is the ratio of the specific heats of the gas. 
Furthermore, if — is a nondimensional parameter, such 
as percentage of chord, designating a surface point on 
the airfoil, then the pressure coefficient, C,, at that 
point is given by 
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The appropriate differential equation for such tran- 
sonic flows is, however, nonlinear and, in general, dif- 
ficult to solve. For the special case of slightly super- 
sonic flow, the problem is, on the other hand, much 
simpler, since the exact solutions associated with this 
problem, the oblique shock and Prandtl-Meyer flow, 
are known. The task is then simply to reduce the 
exact solutions to forms that are appropriate to tran- 
sonic conditions involving the parameter K. This is 
indeed already indicated by von Ka4rman!' himself. 
The purpose of this paper is to complete the calcula- 
tions of von Karman, to render them more systematic, 
and then, finally, to apply these results to calculate the 
rather peculiar aerodynamic characteristics of simple 
thin airfoils in slightly supersonic flow with shocks at- 
tached to the leading edge. 


j (6) 


OBLIQUE SHOCK 


Let 8 be the shock angle—i.e., the angle between the 
free stream and the shock front—and let 6 be the angle 
through which the flow is deflected. The basic oblique 
shock relations’ are then 


r — |] ™ 
| oe M®? sin? g — ~—— (7) 
p® yt1 
and 
tm@-O y-1, 2 1 8 @& 
tan B = y+1 y+1M®? sin? 8 ; 


On the other hand, the pressure coefficient C, for a point 
on the surface is given by 

C, = (2/yM°?*)[(p/p°) — 1] (9) 
In the three equations, Eqs. (7), (8), and (9), the two 


quantities 8 and (p/p°) can be eliminated. The re- 
sultant relation is 


owe. 2 \". 
2\(R +5h3) 7 


oa cs. ~vr 
ali - (F4,)] (10) 


Eq. (10) is the exact oblique shock relation. For the 
transonic conditions, @ is small and is equivalent to 6 
in von Karm4n’s formulas. Thus 


K = (1/2)[(M°? — 1)/(r0)”*] 


tan 6 (1 _ 


(11) 
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From this equation the free stream Mach Number, 
M°, can be expressed in terms of K and 6. In fact, 


1/M°? = 1 — 2K(re@)”* + (12) 


The appropriate form for the pressure coefficient is 
Cy = (0/0) pO(K (13) 


where p"’ is the pressure function associated with the 
flow after the oblique shock. Then Eqs. (12) and (13) 
can be substituted into the general relation, Eq. (10), 
and, by taking only the lowest order of terms, an ex- 
pression for p‘” is obtained: 


(pay?/4) [2K — (p/2)] (14) 


Eq. (14) is a cubic equation for p“. To solve this 


equation, it is convenient to put 
f = 2/p” (15) 


Then Eq. (14) is converted into the normal form of 


cubic equations, 

f? — 2Kf+1 (16) 
According to standard solutions,’® the character of the 
roots of this cubic equation is dependent upon whether 


(1/4) (8/27)K* 


is greater than, equal to, or smaller than zero. In the 
first case, there will be one real root and two complex 
conjugate roots. The complex roots have, of course, 
no physical significance and need not be considered 
further. However, the single root is negative, as can 
be seen from Eq. (16) by making A extremely small. 
Negative f means negative p™. But physically the 
flow is compressed by the shock, and thus p“ must be 
positive. Therefore, no useful solution exists for 
[(8/ /27)K*] > 0 
When K is equal to the critical value A* such that 











TABLE 1 
Pressure Function for W eak Oblique Shock, f:'"), Strong Oblique Shock, 
po), Prandtl-Meyer 2 xpansion, p‘?), and the Ratio of (\f2 — 1)/(M°2 1 
for Various Values of K. [E qs. (19), (20), (29), and ( (36)] 
M2 — 1 M2 _ 
K pil! pr) Me? — |} — pf Moz — :) 
0.00 ° 2 a0 
0.20 l 5.4891 
0.40 1.! 2.9987 
0.60 ; 1 2.2086 
0.80 ‘ 2% a4 1.8333 
0.945 2.5028 2 58 1 1.6682 
0.975 2.1199 3. 2 , 1.3 1.6408 
1.00 1.9993 3.2260 0.0000 1. 1.6193 
1.25 1.5149 t 45 0.3941 1 1.4572 
1.50 1.3051 5.7574 0. 56: 50 Zs 1.3550 
1.75 1.1714 6.8199 0.6 By 1.2859 
2.00 1.0773 7.6295 0.7 0.4 1.2365 
2.50 0.9393 9.9069 0 0.8 1.1717 
3.00 0.8487 11.9204 0 0 1.1318 
4.00 0.7235 15.9490 0.4 0.6 1.0865 
5.00 0.6428 19.6155 0. 0. 1.0623 
6.00 0.5844 23.8720 0.¢ 0.5 1.0475 
8.00 0.5038 32.3625 0.9685 > 0 1.0310 
10.00 0.4496 40.3388 0.9775 0.44 1.0222 
13.00 0.3936 53.0786 0.9849 0 1.0150 
20.00 0.3168 78.4929 0.9921 0 1.0079 
30.00 0.2584 127.7466 0.9957 0 1.0043 
50.00 0.2001 199 .0446 0.9980 0 1.0020 
80.00 0.1581 472.1435 0.9990 0. 1.0010 
100.00 0.1414 oo) 0.9993 0.1420 1.0007 
cy 0 0000 © 1.0000 0.0000 1.0 


1949 


Peat TTT r 
2.5 |+- a oe 
| \ ; | | TA 
“(ly | 
| \ p | Ai 
| ant 





2 








WEAK SHOCK PRESSURE FUNCTION, ef"? 





STRONG SHOCK PRESSURE FUNCTION, 10 


\/K 


Fic. 1. Variation of compression shock pressure function with 
K. [Egqs. (19) and (20).] 
(1/4) — [(8/27)(K*)*] 
or 
(3/2)(1/4'") = 0.945 (17) 


there are two distinct roots: f = and f = — 
The second root is again without tl il siguidicance, 
Corresponding to the first root, the critical value p™ 
is 

pO(K*) = 2°" = 2.502 (18) 


K* is actually then the lowest value for K, or the largest 
value for @ at a given M°, for which the result has physi- 
cal significance. This value of A must then corre- 
spond to the maximum flow deflection or critical condi- 
tion for attachment at the leading edge. This checks 
with the calculation of von Karman. ! 

For larger values of A, the physically meaningful 
values’ of p™ are 


py) = | ¥ K cos e °)| (19) 
pi? = uv \: K cos (r+*)| (20) 


¢ is given by 


and 


where 
cos ¢ = (1/4) (27/2)(1/K®) (21) 


Since K > K*, ¢ is greater than zero but smaller than 
(x/2). Under these conditions, p, > p,™, p,™ 
then corresponds to the weak shock case, and p, 
corresponds to the strong shock case. Following the 
general practice, p,;") will be used for attached leading- 
edge shocks. 

The numerical values of p;" and p.“” are tabulated 
in Table 1, and the data are plotted in Fig. 1. p,"' 
can be expressed as an inverse power series of K. Cal- 
culation shows that 
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1 0.25000 0.11048 
pi? = mak .41422 KA + oa Ki 
0.06250 0.03987 0.02734 
an ar oe 
0.01966 0.01463 + 0.01119 
K? a KY 


0.00873 
Pall + ...¢ (3) 


The first two terms of Eq. (22) agree with von Ka4r- 
man's calculation. He did not consider the additional 


terms. 
Let p* denote the critical pressure corresponding to 


local sonic velocity. Then for isentropic flow, 


1 — afy—t 
y+ “2 (M®? — 1) 
p* 2 2 


p° (y +1)/2 





The corresponding pressure coefficient, C,*, is then, 
according to Eq. (9), 


; 9 y¥-!1 ¥/y-1 \ 
Wo eee a wt 
C 54h t ++1 a » | f 


When \J° = } for transonic flows, the above equation 
simplifies to 


C,* = (M°? — 1)/T (23) 


The flow after the strong shock is always subsonic. 
Sonic velocity after shock can only be associated with 
the weak shock. But it is known that the entropy 
change of a weak shock is negligible, and therefore 
Eq. (23) can be used for the weak shock case. 

In general, as shown by Eq. (13), 


Cy = [(M°? — 1)/2T](pi/K) (24) 


By comparing Eqs. (23) and (24), the critical value of 
(p:)* corresponding to local sonic velocity after the 
shock is 

piP*/K, = 2 (25) 


where A, is the value of K corresponding to (p,“)*. 
By dividing Eq. (14) by K°, 


ey fc nd | PY f ee - 
e731) [2 (4) (26) 


and then substituting Eq. (25) into Eq. (26), the equa- 
tion for K, for sonic velocity after a shock is obtained: 


K, = 1 (27) 


Thus, for supersonic velocity after the leading edge 
Shock, K > 1. 

For airfoil calculations, one will also need the ratio 
(M? — 1)/(M®°? — 1), where M is the Mach Number 
after the shock. This can easily be determined by 
noting that in general for isentropic flow 
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Fic. 2. Variation of (M2 — 1)/(M®? — 1) with (p/K). [Eq. 
(29).] 





| {? + [Cy = D/y + DIG? - JZ 
p° 1+ [(y — D/(y + 1) — 1) 
where both (J/? — 1) and (.V/°? — 1) are small under 
transonic conditions. Then Eq. (9) gives 
C, = [((M®? — 1)/r]}tl — (CA? — 1)/(M®? — 1))} 
28) 

By combining Eqs. (24) and (28), the required ratio is 
given as 

(Af? — 1)/(M°* — 1) = 1 — (1/2)(b/K) (29) 
This formula is, of course, generally true. If the value 
p is taken from Eq. (19), then the ratio will be that 
for the weak oblique shock. The numerical values 
are tabulated in Table 1, and the data are plotted in 
Fig. 2. 


PRANDTL-MEYER EXPANSION 
If g is the velocity and @ is the angle of turn during 
expansion, the Prandtl-Meyer expansion is determined 
by the relation‘ 
dq = qd0/VM? — 1 
However, in general, 
dp = —pq dq 
Thus, by eliminating dg from these two equations, 
dp = —pq’ do/V M? — 1 
or 
(2)-- (2) — 
p° p°*) JM? — 1 
Now since the expansion is adiabatic, 


— | — 1 
a? + —— Gia ti + 2 U2 


Thus by using the isentropic relations, g* can be calcu- 
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lated from the above equation as 


2 =} re 
q’ z= 1 [1 + z — Me? — () | 
. Deg 2 


Furthermore 





dC, 2 


This is the exact equation. It can be simplified by 
noting that for transonic conditions C, is small and 
M® is close to 1. Then, expanding the terms on the 
right-hand side of Eq. (31) into a power series of C,, one 
finds that the significant terms are 

dC, 


if 9 1/2 r : 2 
‘7 = — 2/(M®? — 1) (: — Ve? im - Cy) (32) 


Now, similar to Eq. (24), one can write 


ae, piney (=) (38 
= ——- $54) ,= : B & 
aie if or K 


where p is the pressure function associated with 
the Prandtl-Meyer expansion. Eq. (32) then be- 
comes . 


Diese Oe ("*) 1 
Kk” — —-}=- — (34) 
3 dK\K V2 — (p/K) 


This can easily be integrated. The result is 


[2 — (p/K)]"* = (3/K"*) + C (35) 


where C is the integration constant. To determine the 
constant, note that, when 6 = 0, no change is made in 
the free stream and C, = 0. Then, when K > o, 
p/K — 0. Hence, the constant in Eq. (35) is equal 
to (2)"*. Then Eq. (35) can be solved for p™, and 
the result is 
p = {2K — [(2K)'? + 3]"”} (36) 
Therefore, when 1° = 1, and thus K = 0, the pres- 
sure coefficient is proportional to 6°’, or 
C, = —[(30)"/r'’] (M° = 1) (37) 
This is in agreement with von K4rman’s calculation. ! 
For other values of K, the pressure function, p, and 
the ratio (M? — 1)/(M®? — 1) are tabulated in Table 
| and plotted in Figs. 2 and 3. 


FLAT PLATE AIRFOIL 


The aerodynamic coefficients of an airfoil can be 
obtained by integrating the pressure coefficients over 


yMe? \T 9 
(: Tr" 2 G) y¥—1 
do = M2” 2 ; yM®? _ \-@=D/ 
+} WORE} ——G : 
os yc 
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9 0\7—-1/y 
M?-1= ( —— + M \(e) om &. A 
os p es 


These relations, together with Eq. (9), then convert 
Eq. (30) into 


2 Me? e~i/¢ 
y-1 (: _ 2 G) | 


—— 31) 


+ Me? — 


y+ iy 
y¥-—1 


the surface of the airfoil. If the local surface inclina 
tion to the free stream is 0, the sectional coefficients are 


(1/c) / (Ciiwee CPeazer) OX 


| 





C; = 


A 


( a - ( 1/c) | ( Crrower = CPcgzee) | 
¢c 9 (38) 


(sin @)dx 


Cres = 


(1/c) / (Corower = Capper) 
[C. P. — (c/4) ]dx } 


where c is the airfoil chord and C. P. is the center of 
pressure. 

Consider now a flat plate inclined at an angle, a, to 
the flow such that an oblique shock and a Prandtl- 
Meyer expansion take place at the leading edge of the 
lower and upper surfaces, respectively. By integrating 
Eqs. (38), noting that sin a may be approximated by 
a, and utilizing Eqs. (13) and (33), one has 


ron = (a’*/0'")L(K) ) 
Co = (a’*/T'")L(K) 39) 
—Cmeys = (1/4)C, 


where L(K) = [pi — p™] and K = (M®? — 1)/- 


2(al')”. 
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Fic. 3. Variation of Prandtl-Meyer expansion pressure function 
with K. [Eq. (36).] 
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AIRFOILS IN SLIGHTLY SUPERSONIC FLOW 


Rearranging Eqs. (39), the aerodynamic coefficients 
may be plotted in a simple manner involving only the 


lift function L(K). Thus, 
‘ ‘sam YP r’’ 
L(K) == 2); Cc; = ~ b/3 Ca = —4- 3); Ca 4 (40) 
Qa a Qa 


This has been carried out in Fig. +. Curves of Ci 
versus 1° and C; versus a have been calculated from 
these data and are presented in Figs. 5 and 6. Drag 
and moment coefficient curves can be seen to be similar 
to the lift coefficient curve with only a change in the 
scale of the ordinate. Data obtained by the use of 
Ackeret’s linearized theory have been included with 
the curves for comparison. Two features are immedi- 
ately apparent from the figures: First, the lift versus 
angle of attack curve is not a straight line as predicted 
by the linear theory but has a definite upward curva- 
ture; secondly, the predicted lift is at all times greater 
than Ackeret’s prediction, and this difference increases 
with increasing angle of attack. 


ASYMMETRIC WEDGE AIRFOIL 


As an example of finite thickness airfoils, consider 
now the asymmetric wedge depicted in Fig. 7. Denot- 
ing the vertex angle at the leading edge by 0, Table 2 
lists the C, and K for each region of the airfoil. It 
should be noted that the C, for region (2) given in 
Table 2 is relative to region (1). Account may be taken 
of this by noting that 


\ pe ei p° Pe one p1 Pi _ p® 
(Cp) 2relative eR ee +— 
free stream q q q 
Bs (Cy) 2retative to + (Cp)1 (41 ) 


region (1) 
to the order of magnitude of these calculations. 
Since all angles are small, the coefficients may be 
written in integrated form as 


¢, = (C,)s — (1/2)(Cy)1 — (1/2)(C,)2 

Ce _ a(Cy)s + (1/2) (8 va a)(C,), — (42) 
(1/2)(0 + a)(C,)2 

Can = (1/4) (Ge — (Cre] J 


In order to transform Eqs. (42) into a form similar 
to Eq. (40) for the flat plate, two parameters, K, and 
m are now defined such that 


Ky = (1/2)[(M®? — 1)/(r)”*] (43) 
and 
m = (a/@) (44) 


The expressions for K and C, in terms of these param- 
eters are tabulated in Table 3. Taking into account 
Eq. (41), Eqs. (42) can now be written 
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Lift, drag, and moment similarity curve for flat plate 
airfoil [Eqs. (40).] 
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= mp, (K3) — (1 — m)*"p,(K,) — (2) 7p (Ka) 


= L(k, m) - 
C,(0'*/6" _) = m” *bi) (K3) + (1 —_ m) *b, (Ky) _- (3) (1 + m) p (Ko) (45 
= D(K, m) 9) 


ii Cre A(T 
= M(K, m) 


Aye) = (1/4)[m*"p,\(K3) — (2)°*p(K2)] 





where the right-hand sides of Eqs. (45) are functions 
only of K, and m. These are the similarity relations 
for the asymmetric wedge airfoil. Calculations have 
been carried out and plotted in Figs. 8a, 8b, and 8c, 
for a range of values of 0 Z m Z 1 and 0.632 Z K, Z 
1.90. The chosen upper value of A, was approximated 
by using a free stream with J/° equal to 1.30 as a maxi- 
mum and a wedge with 6 = 3° asa minimum. The 
chosen lower value of K, was determined by allowing a 
small range of values for m in which the present theory 











is applicable. The theoretical value of Ky, —~ ~—i.e., 
either a flat plate or a wedge in a stream of infinite 
M°—is seen to lie along the horizontal axis. Note, 
in Figs. 8a, 8b, and 8c, the limiting line for an attached 
shock (Ks; = 0.945) at the leading edge when approxi- 
mately 0.55 Zm Z 1.0, and note also the limiting 
line for sonic flow (K, = 1.0) in region (1) of the airfoil 
when approximately 0 Z m Z 0.55: The latter is a 
necessary condition for applying the Prandtl-Meyer 
flow in region (2). 
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R = (6°/3/T"/s)p(Ki). 


[Eq. (45).] 
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Determination of Stability Derivatives from 
Flight Data 


Harry Greenberg 
Ames Aeronautical Laboratory, N.A.C.A., Moffett Field, Calif. 


October 13, 1948 


N AN INFORMATIVE ARTICLE in the JOURNAL OF THE AERO- 
NAUTICAL SCIENCES (September, 1947) entitled ‘‘Progress in 

Dynamic Stability and Control Research,’’ by W. F. Milliken, 
Jr., of Cornell Aeronautical Laboratory, it was shown how cer- 
tain dynamic stability parameters of an airplane can be deter- 
mined from oscillation tests made in flight. These consist of 
combinations of the fundamental stability derivatives such as 
Ma, Mg, Ma, etc. In the above-mentioned paper, Milliken and 
his coworkers were able to determine the lift curve slope from the 
flight data but were unable to break down the stability parameters 
b and & into their component parts, the stability derivatives. 
The purpose of this letter is to show that there is a fundamental 
reason why it is impossible in some cases to determine the moment 
derivatives solely frem flight data, no matter how accurately or 
carefully the tests are made. It is possible in these cases to 
determine only certain linear combinations of the derivatives; 
these are sufficient to determine the dynamic characteristics of 
the airplane. 

Consider first the case of longitudinal motion with forward 
speed assumed constant, an assumption that is legitimate at the 
start of a maneuver or for high-frequency oscillation tests. This 
is the case analyzed in some detail in the paper referred to above. 
In this case the equations of motion are: 


aLa = mV(q — a) (1) 
aMa + &Ma + qMy + 6Mg = Iq (2) 


By making simultaneous measurements of @ and g (or of n and 
q, where n is the normal acceleration), it is certainly possible to 
determine La/mV from Eq. (1). It might be thought that, by 
taking a sufficient number of measurements of @ and q, the values 
of Ma/Iy, Ma/Iy, Mq/Iy, and Ms/Iy could be determined from 
Eq. (2) by solving four simultaneous equations for the four 
unknowns. But these equations will not be linearly independent, 
since a linear relation exists between the coefficients of these 
simultaneous equations. This linear relation is contained in Eq. 
(1), which states that the quantities a, g, and @ are linearly de- 
pendent. Hence, only M;/Iy can be determined from Eq. (2). 
As for the other three unknowns (Ma/Iy, Ma/Iy, and M,/Iy), 
only two linearly independent combinations of them can be 
determined. The parameters b and k mentioned in Milliken’s 
report are two such combinations. , 

If the equations are generalized to include the derivatives L; 
and Lg, it is still possible to determine the lift derivatives, but 
again it is not possible to determine the moment derivatives 
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solely from flight measurements involving two degrees of freedom 
for the same reason given above. 

lf the flight measurements are made in such a way as to allow 
the forward speed to vary (low-frequency oscillation tests or long 
duration maneuvers) and include at least two c.g. positions, then 
it is possible to determine the neutral point of the airplane. [If 
the derivative My is zero, then it is possible to determine M, 
and all other moment derivatives. But if the derivative My exists, 
it is again impossible to determine Mg and My individually but 
only a linear combination of them. This can be shown in the same 
manner as above for two degrees of freedom. 

In the case of lateral stability, the situation with regard to the 
determination of the derivatives from flight measurements js 
more favorable. This is because of the assumed absence of a 
yawing moment derivative with respect to rate of change of side- 
slip, analogous to the pitching moment derivative with respect to 
rate of change of angle of attack. Hence, each equation is soly- 
able for its derivatives by the method of simultaneous equations, 
since these equations are linearly independent. 


On Sweat-Cooling 


E. A. Richardson 
Bethlehem Stee! Company, Bethlehem, Pa 
July 31, 1948 


I HAVE GIVEN some study to Preprint No. 137, entitled ‘“Experi- 

mental Study of Cooling by Injection of a Fluid Through a 
Porous Material,’’ by Pol Duwez and H. L. Wheeler, Jr., pre- 
sented at the Third National Flight Propulsion Meeting, March 
19, 1948 (published in the JoURNAL OF THE AERONAUTICAL 
ScrencEs, Vol. 15, No. 9, p. 509, September, 1948). This letter 
will cover not/only this paper, but also certain references as 
“Stability of the Laminar Boundary Layer with Injection of a 
Cool Gas at the Wall,’’ by Lester Lees, Aeronautical Engineering 
Laboratory, Princeton University, May 20, 1948, and a recent 
paper by S. W. Yuan, of the Polytechnic Institute of Brooklyn, 
on the temperature field in the boundary layer. In none of the 
last three papers did the authors appreciate that heat flow ina 
fluid moving with a velocity component parallel to, and opposite 
to, the direction of the flow of such heat is vastly different from 
that in a stationary fluid. In my opinion all such papers are quite 
unrepresentative of the matter in question, however well con- 
ceived otherwise. Reference should be made to the paper “‘Math- 
ematical Attacks cn Forced Convection Problems,” by T. B 
Drew, ‘Transactions of the American Institute of Chemical 
Engineers, Vol. xxvi, 1931, in which the differential equations for 
constant mass velocity of flow are given and the case of the con- 
stant and uniform line source of heat normal to the flow direction 
is derived. A little further study through the derivation of a 
uniform and constant plane source of heat normal to the direction 
of constant mass velocity fluid flow would show that the steady- 
state temperature distribution upstream from such a plane source 
is a long function falling off with extreme rapidity. This tends 
to be true when (cypvx/2k) for the fluid is of the order of, or 
greater than, about 2, where cyp is the specific heat per cubic foot 
at constant pressure, v is the fluid velocity in feet per hour, x is 
the distance from the plane heat source in feet, and & is the con- 
ductivity in B per square foot per hour per 1°F. per foot. Hence, 
with small flow velocities, virtually no heat can travel as far as 
(or the order of) 0.001 ft. upstream. Hence, at attainable flow 
velocities, only radiant heat can pass through a relatively thin 
fluid layer. It should be obvious that a thicker layer, or a faster 
velocity, when the boundary layer varies in velocity normal to 
the source, will yield the same sort of result. 
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Table 3 of Duwez and Wheeler seems incorrect in the quoted 
velocity of fluid flow V of the hot gases. I plotted the heat flow 
qin B per sq.in. per sec. versus V in ft. per sec. and secured sub- 
stantially straight lines for any assigned pair of values of hot 
gas temperature Ty and surface of porous body temperature 7}. 
For Tables 2, 4, and 5 these straight lines fairly consistently passed 
through the origin, but for Table 3 the common point approxi- 
mated g = Oat V = —200 ft. per sec. I then plotted the fluid 
exit velocity v to the gas stream velocity V and again obtained 
substantially straight lines. As before, Tables 2, 4, and 5 gave 
lines through the origin, while Table 3 gave lines through a point 
atv = 0, V = —200 ft. per sec. This matter should be cleared 
up. 

Limiting myself to Tables 2, 4, and 5, I proceeded to derive 
coefficients C’ in the formulas Q = C’V, plotting them versus 
(Tr, — T1;), which last is the difference between the average 
temperature 7z, of the turbulent hot gas core and that, 7}, of the 
hot surface of the porous wall. Three such plots resulted, one for 
each table. It did not appear that they could be reduced to one 
plot through use of Reynolds Number Rg. In each case, the plot 
was effectively a straight line through the origin up to a charac- 
teristic value marking the beginning of a curve convex upwards 
which tended to exhibit a maximum. Lack of range in the data 
resulted in but cne curve extending far enough to show a maxi- 
mum followed by a rapid falling away therefrom, this curve being 
for Tz = 1,900°F., which showed a maximum at approximately 

Tr — 71) = 1,400°F. 

In similar fashion, I then plotted coefficients C” in formulas 
of the type v = C”V against (Tz — 7) for each value of Ty. 
The resulting plots in each case comprised two intersecting 
straight-line sects. The point of intersection of these sects corre- 
sponds approximately to the point of departure of C’ from the 
straight line in the preceeding set of plots. 


It will be seen, I think, that any attempt to study the ratio of 
T — T:)/(T: — Te) (where T is the hot gas temperature and 7) 
that of the entering fluid) against Q/W, the ratio of permeating 
fluid flow per square inch per second to that of the hot gases, as 
in Fig. 10 where a narrow and straight band of points passing 
through the origin is found from 0 to 0.014, can be misleading. 
Similarly, the plotting, following Rannie, of (T — 7))/(T — To) 
which should yield curves asymptotic to 1.0, as Q/W increases, 
is misleading, since even the plot of experimental points in Fig. 
11 show a definite trend to a straight line (or lines) above Q/W = 
0.011, which would reach 1.0 by, or before, Q/W = 0.030. This 
isextremely important if true. The plotting of (T — 7)/(T, — 
Ts) versus v/U, as in Fig. 12, is not so misleading, since the de 
partures from the mean (or extrapolated) curve are noteworthy 
above v/U = 0.0047. It is of interest that similar points tend to 
lie on straight lines above this point. (U differs from V in that 
the former is the flow average velocity for the stream while the 
latter is the maximum velocity measured at the core center.) 
This curve may well be extremely misleading, however, since the 
form of plotting used by the present writer indicates the easy 
attainment of much higher ratios of (T — 7\)/(T; — T >) with 
moderate increases of v/U and would look like the stress curve 
in column buckling. (See Timoshenko, Theory of Elastic Sta- 
bility, Ist Ed., p. 39, Fig. 18; McGraw Hill Book Company, Inc., 
New York, 1936). 

With reference to Table 1, in which nitrogen and hydrogen 
cooling are compared, the ratio of the flow of the two fluids for 
equivalent results is that demanded by the theory of heat diffusion 





in that the fluid exit velocity ratio must be that of V bee; /RiCo 
where k is the conductivity, c is the heat capacity of the fluid, 
and subscript 2 might refer to hydrogen and 1 to nitrogen. As 
4 is substantially equal to c2, for these fluids, the ratio is that of 
the conductivities. The calculated ratio is closely the measured 
ratio. My original calculation was too low because of the use of 
old data; more modern data give a better check. 


FORUM 63 


My own experience in the art of cooling high-temperature 
furnaces indicates that it is not necessary to use a separate fluid 
for a coolant. The desired wall protection may be secured by 
passing the whole air mass for combustion through the furnace 
walls, if desired, since it must be obvious, from the same art, that 
vortex-type combustion furnaces with the highest rates of heat 
generation may be used because this same boundary layer keeps 
the burning gases from licking the wall. 

Although heat flow by radiation from the gases was small in 
the experimental apparatus (of the order of 0.001B per sq.in. per 
sec.), the same will not necessarily be true in a larger apparatus, 
particularly when combustion is occurring therein. 


- 





Induction Effects in Aeroelasticity 


John W. Miles 
Department of Engineering, University of California at Los Angeles 


September 7, 1948 


ILDEBRAND AND REISSNER have investigated the effects of 

induction on the torsional divergence speed of an unswept 
wing and concluded that strip theory, using the two-dimensional 
lift curve slope, led to errors of from 15-40 per cent.!. The writer 
recently had occasion to ifivestigate this problem for Northrop 
Aircraft, Inc.,2 and found that the error introduced by strip 
theory for elliptic or tapered wings of continuous stiffness dis- 
tribution is made small simply by introducing the mean lift 
curve slope for the rigid wing (i.e., corrected for induction effects) 
rather than the two-dimensional slope. 

Analysis of the integrodifferential equation governing the 
torsional deformation of the aerodynamically loaded wing leads 
to the introduction of the modulus 6 as the ratio of the indicated 
air speed to the divergence speed (i.e., @ = 1 at divergence), 
and it is found that 


1 
o = ais] ff e(s)c#(s)0(s)c°(s)ds | x 
0 
"4 fe) ]2 ) 
1 I cro] 252 as} (1) 
/J0 ) 


where g is the dynamic pressure, 0 is the span, e is the fraction of 
the chord of the elastic axis aft of the line of aerodynamic centers 
(eccentricity), c is the chord, @ is the elastic twist, @° is the lift 
coefficient distribution due to 6, GJ is the torsional stiffness 
distribution, and s is a dimensionless coordinate running from 0 
at the root to 1 at the tip. Morecver, the equation is self- 
adjoint, & being the eigenvalue, and it therefore follows (the Ray- 
leigh principle) that ® is an absolute minimum with respect to 
first order variations of 6(s) about the true solution to the integro- 
differential equation. 

Using the Prandtl lifting line theory to determine cy” from @, 
using variational (or Rayleigh) theory to solve Eq. (1), and as- 
suming that GJ(s) was proportional to the fourth power of the 
chord, the writer found? that @ was given by 

& = (3/8) [(dez/da)eS%g/GI(0)}'/* (2) 
with a maximum error of 3 per cent for taper ratios between 1.6:1 
(dcz/da) is the lift curve slope for the rigid wing, and e 
is assumed constant. For asymmetrical loading (dcz/da) should 
be evaluated for half the true aspect ratio. Thus, it is well to 
introduce @; and &q for symmetrical and asymmetrical loading. 


and a. 


1 Hildebrand, F. B., and Reissner, Eric, The Influence of the Acrodynami 
Span Effect on the Torsional Divergence Velocity and on the Shape of the 
Corresponding Deflection Mode, N.A.C.A. T.N. No. 926, 1944. 

2 Miles, John W., Aeroelastic Properties of Unswept Wines, Lifting Line 
Theory, Northrop Report A-101, 1947 
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The modulus @ is the dominant parameter in all aeroelastic 
calculations not involving control surface loading. Thus, it is 
found? that the lift coefficient (given root angle of attack) and 
damping in roll coefficients for a wing of two to one taper and as- 
pect ratio 6 (and, within a few per cent, for all unswept wings) 
are given by 

(cr/ex®) = 1 + 0.87 [a52/(1 — &s2)] (3) 

(ctp/cip?) = 1 + 0.98[Gg?/(1 — ca?)] (4) 
where the superscript 0 denotes the value of the coefficient for the 
rigid wing. 

In the case of control surface loading the reversal speed is more 
significant than the divergence speed, and it is therefore expedient 
to introduce the parameter @p as the ratio of the indicated air 
speed to the reversal speed (@g = 1 at reversal). In this case it is 
found that the results are relatively insensitive to induction 
effects if the two-dimensional lift curve slope is used. The ratio 
(@p/@), as given by Eqs. (1) or (2) [with (dcz/da) evaluated as 
the two dimensional lift curve slope], is plotted against the 
parameter, I's in the accompanying Fig. 1 for a 30 per cent span 
aileron on a two to one tapered wing of aspect ratio 6, where 


1949 


rs =1+ (Cm, /eCts) (5) 
ms and cls being the section characteristics of the control surface. 
These results are not too sensitive to variations of parameters 
other than Is. Using N.A.C.A. test data,* I's is plotted versus 
the control surface chord ratio for various products of the two 
dimensional lift curve slope (m) and the eccentricity (e) in Fig. 2, 
If Sp is evaluated according to the last paragraph, it is found 
that 
(p/p®) = 1 — ap? (6) 
to within a few per cent for all practical configurations. This 
parabolic fall-off of rolling effectiveness with indicated air speed 
seems to be a rather well-known empirical conclusion. 

More detailed results, including the calculation of dynamic 
stability derivatives, are given in reference 2. It should be spe- 
cifically remarked that theresults are not applicable to swept wings, 
where bending will generally be more important than torsion. 

+ Ames, M. B., and Sears, R. I., Determination of Control Surface Char- 
acteristics from N.A.C.A. Plain Flap and Tab Data, N.A.C.A. T. R. No. 721, 
1941. 





Airfoils in Slightly Supersonic Flow 


(Continued from page 61) 


To illustrate the nature of the coefficients in the 
usual form, calculations have been carried out for the 
case of 9 = 3° and are plotted in Figs. 9 through 12. 


Fig. 9 shows that, at small angles of attack giving 
positive lift, the lift coefficient actually increases with 
increase in Mach Number. This is a behavior not 
predicted by Ackeret’s linear theory, which requires a 
decrease in lift coefficient with Mach Number according 
to the factor 1/(M°? — 1)'”. 
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